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Abstract 

A  field-deployable  hyperspectral  chromotomographic  imager  has  been  developed 
and  tested  as  a  risk-reduction  prototype  to  assist  design  of  a  space-based  system. 
The  instrument  uses  a  high-speed  video  camera  looking  through  a  rotating  direct- 
vision  prism  to  simultaneously  observe  the  full  field  of  view  in  all  visible  wavelength 
channels.  This  enables  hyperspectral  imaging  of  transient  events  at  high  temporal 
resolution.  The  chromotomographic  process  multiplexes  the  spectral  and  spatial  so 
an  advanced  reconstruction  algorithm  is  required  to  separate  the  spectral  channels. 

A  physics-based  model  of  the  instrument  was  developed  to  assist  in  future  trade- 
space  choices  for  design  of  the  spaced-based  system.  Additionally,  the  model  is  used 
for  the  development  and  assessment  of  the  filtered  backprojection  reconstruction  algo¬ 
rithm.  The  model  is  used  to  simulate  various  scenarios  and  gauge  the  reconstruction 
algorithm’s  performance.  An  assessment  on  how  the  reconstruction  algorithm  per¬ 
forms  under  chromatic  and  monochromatic  aberrations  is  included. 

Laboratory  experiments  from  the  field-deployable  instrument  were  collected,  and 
the  results  are  compared  to  physics-based  model  predictions.  Prior  to  reconstructing 
scenes  collected  by  the  field  instrument,  a  transverse  offset  resulting  from  an  imper¬ 
fection  of  the  direct-vision  prism  was  characterized  and  corrected.  Results  from  the 
simulated  and  experimental  data  show  that  the  instrument  and  algorithm  are  capable 
of  detecting  spectral  and  spatial  information  of  complex  scenes. 
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CHARACTERIZATION  OF  A  HYPERSPECTRAL  CHROMOTOMOGRAPHIC 


IMAGING  GROUND  SYSTEM 

1.  Introduction 

This  thesis  presents  a  physics-based  model  and  experimental  validation  of  the 
ground-based  chromotomographic  experiment  (GCTEx)  developed  by  the  Air  Force 
Institute  of  Technology  (AFIT).  The  optical  model  is  constructed  using  Zernax  and 
MATLAB.  The  model  provides  the  mechanism  required  to  analyze  future  trade-space 
engineering  choices  for  the  design  of  the  space-based  system. 

Previous  work  done  at  AFIT  resulted  in  models  [9,  28]  and  reconstruction  algo¬ 
rithms  [16,  18]  using  ideal  optics.  However,  it  is  apparent  that  the  GCTEx  instrument 
is  not  ideal  and  optical  aberrations  coupled  with  systemic  errors  are  present  [35-37]. 
These  systemic  errors  and  aberrations  limit  the  spectral  and  spatial  performance  of 
the  GCTEx. 

There  are  three  primary  objectives  of  this  research:  development  of  a  physics- 
based  model  that  allows  for  controlled  amounts  of  systemic  error  and  aberrations, 
the  evaluation  of  how  systemic  errors  and  aberrations  impact  the  reconstruction  al¬ 
gorithm,  and  experimental  validation  of  the  physics-based  model. 

1.1  Motivation 

This  thesis  focuses  on  one  field  of  remote  sensing  known  as  hyperspectral  imager 
(HSI).  Hyperspectral  imaging  detects  both  the  spatial  and  spectral  content  from  an 
observed  scene.  Conventional  HSI  devices,  such  as  scanned-slit  or  optical-filtered, 
are  well-suited  to  observe  static  or  slowly  varying  scenes  and  are  commonly  used 


1 


in  air  and  space  surveillance.  However,  these  HSI  devices  have  a  limited  ability  to 
observe  scenes  that  vary  quickly.  A  capability  gap  throughout  HSI  devices  exists  for 
capturing  all  the  spatial  and  spectral  information  of  transient  events  (phenomenology 
varying  more  than  10  Hz).  One  application  of  collecting  all  aspects  of  transient  events, 
such  as  bomb  detonations,  allows  quick  battle  damage  assessment  or  detection  for 
secondary  detonations.  Using  a  direct-vision-prism  (DVP)  and  a  high-speed  camera, 
hyperspectral  chromotomographic  imagers  simultaneously  observe  the  full  field-of- 
view  (FOV)  in  all  spectral  channels  necessary  to  capture  these  transient  events. 

1.2  Hyperspectral  Remote  Sensing 

Hyperspectral  remote  sensing  combines  two  sensing  modalities:  imaging  and  spec¬ 
trometry.  An  imaging  system  captures  a  picture  of  a  remote  scene  related  to  the 
spatial  distribution  of  the  power  of  reflected  and/or  emitted  electromagnetic  radi¬ 
ation  integrated  over  some  spectral  band.  Spectrometry  measures  the  variation  in 
power  with  the  wavelength  or  frequency  of  light,  which  captures  information  related 
to  the  chemical  composition  of  the  materials  measured.  The  instrumentation  used 
to  capture  such  spectral  information  is  called  a  spectrometer.  Hyperspectral  imaging 
sensors  are  able  to  simultaneously  capture  both  the  spatial  and  spectral  content  of 
remote  scenes. 

Hyperspectral  imagers  collect  radiation  from  the  electro-optic  (EO) /infrared  (IR) 
spectral  regions,  defined  as  the  a  portion  of  the  electromagnetic  spectrum  that  nom¬ 
inally  ranges  from  0.4  to  14  micron  wavelength  (20  to  750  THz  frequency).  The  elec¬ 
tromagnetic  radiation  is  often  segmented  into  five  spectral  regions:visible  (VIS)  from 
0.4  to  0.7  pm,  near  infrared  (NIR)  from  0.7  to  1.1  pm,  shortwave  infrared  (SWIR) 
from  1.1  to  3.0  pm,  midwave  infrared  (MWIR)  from  3  to  5  pm,  and  longwave  in¬ 
frared  (LWIR)  from  5  to  14  pm.  The  spectral  region  in  which  hyperspectral  imagers 
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operate  is  dependent  on  the  object  of  interest. 

For  the  purposes  of  this  thesis,  all  objects  will  be  considered  passively,  vice  actively, 
illuminated.  Passive  objects  use  two  primary  sources  of  radiation:  reflected  light 
and/or  self  emission  from  the  objects  in  the  scene.  Bomb  detonations  or  muzzle  flashes 
are  examples  of  passive  sources  that  self-emit.  Reflected  sunlight  is  the  dominant 
radiation  source  in  the  VIS,  NIR,  and  SWIR,  while  thermal  emission  dominates  in 
the  LWIR.  Key  components  of  the  HSI  are  selected  based  on  what  spectral  region 
the  object  is  best  observed  at.  In  general,  it  more  difficult  to  capture  targets  in  the 
LWIR  region  vice  the  VIS.  The  primary  cause  is  detectors,  or  focal  plane  arrays 
(FPAs),  made  for  the  LWIR  are  sensitive  to  objects  that  have  temperatures  near 
room  temperature.  The  FPA,  imaging  optics,  and  support  equipment  must  be  cooled 
such  that  it  does  not  generate  signal  to  overwhelm  signal  produced  by  the  scene. 

Hyperspectral  imagers  produce  data  commonly  referred  to  as  a  hypercube,  which, 
as  depicted  in  Figure  1.1(a),  is  a  three-dimensional  data  set  composed  of  layers  of 
grayscale  images. 


Figure  1.1.  (a)  Hypercube  taken  from  AVIRIS  platform  over  Moffet  Field,  CA.  (b) 
Example  of  the  spectrum  obtained  from  a  single  pixel  of  the  hypercube.  [26] 
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As  shown  in  Figure  1.1(b),  each  pixel  of  the  hypercube  contains  a  finely  sampled 
spectrum  whose  features  are  related  to  the  materials  contained  within  it.  Each  layer 
is  a  spectral  band  is  referred  to  as  a  bin.  The  ability  to  capture  a  continual  spectrum 
as  opposed  to  more  discrete  set  of  often  discontinuous  spectral  bands  distinguishes  a 
hyperspectral  imager  from  a  multispectral  imager. 

1.3  Research  Goals 

This  research  complements  previous  work  undertaken  by  AF1T  which  has  the 
end  goal  of  fielding  a  hyperspectral  chromotomographic  imager  (CT1)  aboard  the 
International  Space  Station  (ISS).  This  thesis  focuses  on  developing,  evaluating,  and 
validating  a  physics-based  model  suitable  for  future  trade-space  choices  in  the  design 
of  the  space-based  system.  Furthermore,  this  research  demonstrates  how  optical 
aberrations  and  systemic  errors  impact  the  spatial  and  spectral  resolution.  Specific 
goals  addressed  in  this  thesis  include: 

•  Further  progress  the  development  of  the  ground-based  CTEx  instrument 

—  Fix  the  electronics  system  to  accurately  log  prism  angles  with  camera  data 

—  Automate  the  reconstruction  process,  no  manual  corrections  to  produce 
hypercubes 

•  Develop  model  for  trade  studies  for  space  system 

•  Investigate  the  effect  of  aberrations  impact  the  system  performance 

•  Determine  the  operating  modes/limits  of  the  system  and  the  reconstruction 
algorithm 

•  Collect  laboratory  and  held  data  to  demonstrate  the  utility  of  a  CTI 

1.4  Organization 

Chapter  2  provides  an  overview  of  the  optical  concepts  and  that  are  used  through¬ 
out  the  thesis.  Furthermore,  Chapter  2  also  tackles  the  analytical  expressions  that 
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govern  the  device  and  an  overview  of  the  reconstruction  algorithm.  Chapter  2  con¬ 
cludes  with  a  detailed  overview  of  AFIT’s  GCTEx.  Chapter  3  is  a  development  of 
the  physics-based  models  made  from  the  combination  of  Zemax  and  MATLAB.  The 
validation  of  the  models  and  the  source  of  the  transverse  offset  is  also  explained.  Chap¬ 
ter  4  is  an  assessment  of  the  reconstruction  algorithm  from  a  variety  of  simulations. 
Laboratory  data  collected  from  the  GCTEx  is  also  presented.  Chapter  5  concludes 
this  thesis  summarizing  key  results  and  recommendations  for  future  research. 
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2.  Background 


This  chapter  provides  an  overview  of  the  concepts  used  to  construct  the  simulation 
and  reconstruction  algorithm.  The  beginning  sections  cover  the  foundation  of  linear 
system  theory,  basic  principles  of  Fourier  optics,  and  how  aberrations  impact  the 
optical  system.  The  following  section  is  a  brief  overview  of  conventional  hyperspectral 
imagers  and  how  a  trade  study  of  how  they  are  capable  of  measuring  transient  events 
compared  to  CTIs.  The  principles  of  the  hyperspectral  chromotomographic  imager 
operation,  development  reconstruction  algorithm,  and  its  performance  limits  are  then 
covered.  Finally,  a  detailed  overview  the  AFIT  GCTEx  instrument  is  given. 

2.1  Linear  Systems  Theory  of  Optical  Systems 

These  two  sections  briefly  cover  the  basic  principles  of  optical  linear  system  theory 
and  Fourier  optics.  See  the  texts  by  Holst  and  Lomheim  [23]  and  Goodman  [15]  for 
a  thorough  treatment.  Linear  system  theory  was  developed  for  electronic  circuitry 
and  has  been  extended  to  optical,  electro-optical,  and  mechanical  systems  [23].  For 
a  system  to  be  linear,  it  must  meet  the  principles  of  superposition  and  homogeneity. 
Optical  linear  systems  make  use  of  the  property  known  as  shift  invariance  to  calculate 
how  images  are  formed  from  objects.  Shift  invariance  states  that  a  shift  of  the  input 
causes  a  corresponding  shift  in  the  output.  For  electronic  circuits,  the  shift  is  in 
time,  whereas  for  optical  systems,  the  shift  is  in  space.  For  example,  with  an  imaging 
system,  as  a  target  moves  from  the  top  of  the  field  of  view  to  the  bottom,  the  image 
also  moves  from  the  top  to  the  bottom. 

Goodman  showed  that  when  an  optical  system  produces  an  image  using  incoherent 
light  (the  phase  of  the  source  is  a  random  variable),  then  the  function  which  describes 
the  intensity  in  the  image  plane  produced  by  a  impulse  in  the  object  plane  is  called 
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the  point  spread  function  (PSF)  [15].  The  object  is  represented  as  a  sum  of  impulses 
within  its  boundaries.  The  object,  o(x,y),  is  modeled  as  a  series  of  weighted  Dirac 
delta  functions.  For  small  increments,  the  response  to  these  impulses,  h(x,y ),  form 
the  image  intensity  distribution,  i(x,y)  given  by  the  convolution  integral 

/oo  roo 

/  o(x',y')h{5(x  -  x')5(y -y')}dx'dy'.  (2.1) 

-OO  J  —  OO 

Equation  (2.1)  is  symbolically  represented  by  the  two-dimensional  convolution  oper¬ 
ator, 

i(x,y)  =  o{x,y)®h{x,y),  (2.2) 

and  by  transforming  to  the  Fourier  domain  the  convolution  becomes  a  multiplication 

I(u,v)  —  0(u,v)H(u,v).  (2.3) 

The  function  H(u,v )  is  called  the  Transfer  Function  and  in  the  case  of  optical 
systems,  it  is  known  as  the  Optical  Transfer  Function  (OTF).  The  OTF  is  a  complex 
valued  function,  composed  of  real  and  imaginary  parts: 

OTF{u,v )  =  \H(u,v)\ej^u’v\  (2.4) 

The  real  portion  is  the  modulus  of  H(u,  v )  and  is  called  the  Modulation  Transfer 
Function  (MTF).  The  MTF  is  the  ratio  of  output  modulation  to  input  modulation 
normalized  to  unity  at  zero  frequency.  The  phase,  v )  is  termed  Phase  Transfer 
Function  (PTF).  The  impact  of  the  PTF  is  significant  to  the  response  of  the  op¬ 
tical  system;  however,  the  MTF  is  more  effective  in  describing  the  optical  system’s 
performance. 
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2.1.1  Modulation  Transfer  Function. 


Modulation  Transfer  Functions  provide  an  objective  mechanism  for  describing  the 
optical  system’s  ability  to  resolve  spatial  content.  Additionally,  individual  subsys¬ 
tems,  such  as  lenses  and  mirrors,  have  their  own  MTF  which  can  be  cascaded  with  all 
other  the  system  components  to  yield  the  system  MTF.  Figure  2.1  is  a  normalized, 
diffraction-limited  MTF  of  a  symmetrical  optical  system  with  a  circular  aperture. 


Figure  2.1.  MTF  of  a  diffraction-limited  system  with  frequency  normalized  to  the  cutoff 
frequency,  assuming  a  circular  aperture. 


The  horizontal  axis  indicates  the  spatial  frequency,  typically  specified  in  lines/mm, 
for  which  the  optical  system  can  resolve.  The  vertical  axis  is  the  amount  of  modula¬ 
tion,  or  contrast,  defined  as 


M  odulation 


Imax  Imin 
Imax  T  Imin 


(2.5) 


The  variables  Imax  and  Imin  are  the  maximum  and  minimum  signal  levels,  respectively. 


Imax  is  the  maximum  intensity  produced  by  an  image  (white)  and  Innn  is  the  minimum 
intensity  (black).  As  the  frequency  increases,  the  amount  of  modulation  decreases, 
assuming  a  clear  aperture  and  no  apodization  is  applied.  The  highest  input  frequency 
that  can  be  resolved  is  defined  as  the  cutoff  frequency,  fc.  For  a  circular  aperture  and 
incoherent  imaging,  the  cutoff  frequency  is 


fc  = 


D 


1 


(2.6) 


A  Uft  A  F/#’ 

where  feff  is  the  optical  system  effective  focal  length  and  D  is  the  entrance  pupil 
diameter.  The  F-number,  F1/^,  is  the  ratio  of  the  effective  focal  length  over  the 
entrance  pupil 


F/#  = 


feff 

D  ' 


(2.7) 


Because  the  cutoff  frequency  is  wavelength  dependent,  polychromatic  MTFs  are  cal¬ 
culated  individually  for  all  the  wavelengths  of  interest. 

A  diffraction-limited  system  is  the  theoretical  best  performance  for  an  optical 
system  to  achieve.  If  a  point  source  was  imaged  by  a  diffraction-limited  system  with 
a  circular  aperture,  an  Airy  disk  pattern  would  form  on  the  image  plane.  The  Rayleigh 
Criterion  [22]  states  that  this  in  order  to  resolve  two  points,  the  minimum  distance, 
r,  between  the  peaks  is  given  by 


r  =  1.22  =  1.22AF1/#.  (2.8) 

Aberrations  within  an  optical  system  will  degrade,  and  thereby  lower  the  MTF 
response  below  the  diffraction  limit.  Aberrations  may  have  a  significant  impact  on 
the  performance  of  the  optical  system. 
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2.2  Aberrations 


Aberrations  are  the  result  of  optical  components  having  non-ideal  characteristics 
from  those  set  by  Gaussian  optics.  Gaussian  optics  make  use  of  the  paraxial  ap¬ 
proximation  to  assume  that  all  light  rays  from  one  point  of  the  object  through  the 
optical  system  converge  to,  or  diverge  from,  a  single  point.  Aberrations  fall  into  two 
primary  categories:  chromatic  and  monochromatic.  Monochromatic  aberrations  are 
often  associated  with  the  geometry  of  the  optical  component  resulting  in  rays  that 
do  not  converge  to  a  single  point.  Seidel  aberrations  are  often  used  to  describe  these 
monochromatic  aberrations.  Chromatic  aberrations  are  the  result  of  the  wavelength 
dependence  in  the  index  of  refraction.  Chromatic  aberrations  primarily  exist  in  only 
lens-based  optical  systems.  The  following  sections  provide  an  overview  of  aberrations, 
see  the  text  by  Hecht  [22]  for  in-depth  coverage. 

2.2.1  Chromatic  Aberrations. 

The  velocity  of  light  changes  when  it  passes  through  different  mediums.  Short 
wavelengths  travel  slower  in  glass  than  long  wavelengths  causing  the  colors  to  dis¬ 
perse.  Chromatic  aberrations  are  the  result  of  this  dispersion.  In  some  applications, 
this  dispersion  is  desirable  such  as  prisms.  However,  for  most  image-forming  optical 
systems,  chromatic  aberrations  cause  the  focus  to  vary  as  a  function  of  wavelength. 
This  usually  results  in  the  optical  system  to  be  in  focus  for  a  narrow  band  of  wave¬ 
lengths  leaving  the  other  wavelengths  to  be  defocused. 

Achromatic  lenses  are  a  method  of  correcting  for  chromatic  aberrations.  By  com¬ 
bining  sets  of  low  dispersion  and  high  dispersion  lenses,  the  effects  of  the  chromatic 
aberrations  can  be  reduced.  However,  it  becomes  exceedingly  difficult  to  design  an 
achromatic  lens  systems  over  a  large  spectral  range.  Commonly  mirrors  are  used  to 
avoid  the  detrimental  effects  of  chromatic  aberrations. 
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2.2.2  Seidel  Aberrations. 


Seidel  aberrations  are  the  result  of  incorporating  an  additional  term  from  the 
Taylor  series  expansion  made  in  the  paraxial  approximation.  There  are  five  Seidel 
aberrations:  spherical,  coma,  astigmatism,  held  curvature,  and  distortion.  Defocus 
and  tilt  are  first-order  aberrations  often  included  in  the  Seidel  aberrations. 

Seidel  aberrations  are  useful  metrics  for  designing  optical  systems  to  gauge  system 
performance.  The  Seidel  aberrations  are  a  function  of  object’s  location  in  terms  of 
exit  pupil  coordinates,  polarmetrically  described  by  p  and  6 ,  along  with  a  field-height 
dependence,  h.  The  Seidel  aberrations  are  summarized  in  Table  2.1. 


Table  2.1.  First-order  and  Seidel  Aberrations. 


Description 


Focus 

Tilt 


Spherical  Aberration 
Coma 

Astigmatism 
Field  Curvature 
Distortion 


Term 

iW 

Wmhp  cos(0) 
W040P4 

W131  hp3  cos(6>) 
W222h2p2  cos(d)2 
W220h2p2 
W3Uh2p2  cos (0)2 


The  coefficients  in  front  of  each  term  describe  the  severity  of  the  aberration  on  the 
system.  The  weight  of  the  Seidel  coefficient  is  based  upon  the  optical  path  difference 
caused  by  the  aberration  in  terms  of  the  wavelength  of  interest.  It  is  often  convenient 
to  measure  and  describe  aberrations  in  terms  of  a  wavefronts.  A  method  for  describing 
the  shape  of  the  wavefront  is  known  as  the  Zernike  Polynomials  discussed  in  the  next 
section. 
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2. 2. 2.1  Zernike  Polynomials. 


The  Zernike  polynomials  are  a  set  of  orthogonal  functions  that  are  useful  for 
representing  wavefront  aberrations.  They  describe  most  real-world  optical  surfaces 
or  wavefronts  to  be  represented  as  a  short  string  of  coefficients.  Interferometers  often 
output  their  measurements  in  terms  of  Zernike  polynomials.  Additionally,  they  have 
convenient  mathematical  properties  [45]: 

•  Continuous  and  orthogonal  when  used  over  a  circular-shaped  domain 

•  Complete  basis  set,  any  well-behaved  function  is  described  by  a  linear  combi¬ 
nation  of  Zernike  polynomials 


•  Rotationally  symmetric 

•  Computationally  efficient  being  polynomial  in  form 


Zernike  polynomials  arise  in  the  expansion  of  a  wavefront  function  for  optical 
systems  with  circular  pupils.  The  odd  and  even,  respectively,  Zernike  polynomials 
are  given  by  [42] 

Z~m(p,  6)  sin 

=  K(p)  (md),  (2.9) 

Zn(p,d)  COS 

where  9  is  the  azimuthal  angle  defined  within  0  <  6  <  2tt  and  p  is  the  radial  distance 
bounded  between  0  <  p  <  1.  The  radial  function  K'n(p)  is  defined  for  n  and  m 
integers  with  n  >  m  >  0  by 


K(p) 


y-»(n-m)/2  (-phn-Q! _  n-2 1 

Z^l=0  ;!(|(n+m)  —  i)!(|(n— m)— 1)1  ' 


for  n  —  m  even 


0 


for  n  —  m  odd 


(2.10) 


Equations  (2.9)  and  (2.10)  are  evaluated  to  produce  the  commonly  used  FRINGE 
set  of  Zernikes  that  is  composed  of  37  terms.  The  FRINGE  polynomials  a  subset 
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of  the  standard  Zernike  polynomials  developed  by  Born  and  Wolf  [3]  that  have  been 
arranged  in  a  different  order.  Only  the  first  nine  terms  are  listed  in  Table  2.2  and 
like  the  Seidel  aberrations  they  are  represented  polarmetrically  of  the  pupil. 


Table  2.2.  First  nine  Zernike  terms  used  to  represent  the  first  and  third-order  aberra¬ 
tions.  [43] 


Term  # 

n 

m 

Polynomial 

Description 

0 

0 

0 

1 

Piston 

1 

1 

1 

p cos(d) 

X-tilt 

2 

1 

-1 

psin(d) 

Y-tilt 

3 

2 

0 

2 p2  -  1 

Focus 

4 

2 

2 

p2  cos(2d) 

Astigmatism  @  0°  &  Focus 

5 

2 

-2 

p2  sin(2d) 

Astigmatism  @  45°  &  Focus 

6 

3 

1 

p(3p2  —  2)  cos(0) 

Coma  &  X-tilt 

7 

3 

-1 

p(3p2  —  2)  sin(0) 

Coma  &  Y-tilt 

8 

4 

0 

6p4  —  6p2  +  1 

Spherical  &  Focus 

In  Table  2.2,  the  first  3  terms  represent  the  first-order  terms  associated  with 
Gaussian  optics.  Term  0  is  a  constant  or  piston  term,  while  terms  1  and  2  are  tilt 
terms  and  term  3  represents  focus.  The  third-order  aberrations  are  represented  by 
terms  4  through  8.  Terms  4  and  5  are  astigmatism  plus  defocus  and  terms  6  and 
7  represent  coma  plus  tilt,  while  term  8  represents  third-order  spherical  with  focus. 
Terms  9  or  higher  represent  the  higher-order  aberrations. 

Wyant  showed  that  by  using  a  combinations  of  the  first  nine  Zernike  terms,  the 
Seidel  aberrations  can  be  represented  and  are  shown  in  Table  2.3  [43].  However, 
the  Zernike  terms  have  no  dependence  on  held  height  so  the  representations  are  not 
true  Seidel  terms.  Thus  tilt  and  distortion  are  indistinguishable  from  one  other. 
Additionally,  held  curvature  appears  to  be  focus,  or  vice  versa. 

The  Zernike  aberrations  are  used  to  add  controlled  amounts  of  aberrations  into 
the  simulation.  Since  Seidel  aberrations  are  often  used  to  evaluate  optical  system 
performance,  the  Zernike  combinations  listed  in  Table  2.3  are  used. 
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Table  2.3.  Zernike  combinations  equivalent  to  Seidel  Terms.  The  Zernike  combinations 


[escribe  the  magnitude  and  the  phase  terms  have  been  excluded. 

Seidel  Coefficient 

Zernike  Combination 

Description 

Hii 

y/{Z1  -  2Ze)2  +  (Z2  -  2 Z7)2 

Tilt 

w20 

2Z3  -  6Z8  ±  y/Zj  +  Z$ 

sign  chosen  to  minimize  absolute  value  of  magnitude 

Focus 

w22 

±2  y/Zl  +  Z\ 

sign  opposite  that  chosen  in  focus  term 

Astigmatism 

Hb, 

3  y/Z2  +  Z? 

Coma 

w40 

6  Z8 

Spherical 

2.3  Hyperspectral  Imagers 

Hyperspectral  imagers  are  broken  into  two  categories:  dispersive  and  interfero¬ 
metric.  Each  type  has  their  own  advantages  which  is  dependent  on  the  application. 
A  brief  overview  on  conventional  HSI  instruments  and  a  trade  study  on  how  capable 
each  are  on  capturing  transient  events  is  now  given. 


2.3.1  Fourier  Transform  Spectrometers. 

A  common  interferometric  HSI  is  the  Michelson  interferometer  based  Fourier 
Transform  Spectrometer  (FTS).  One  path  of  the  Michelson  interferometer  is  varied  in 
length  with  respect  to  another  using  a  beamsplitter  and  a  set  of  mirrors.  Figure  2.2 
illustrates  the  hypercube  generated  from  the  FTS.  As  the  mirror  sweeps,  thereby 
changing  the  optical  path  length,  the  resulting  images  of  the  changing  interference 
patterns  are  then  correlated  with  the  Fourier  transform  of  the  spectrum  of  the  scene. 

Fourier  Transform  spectrometers  offer  high  spectral  resolutions  but  require  time 
to  sweep  the  mirror.  Additionally,  FTSs  require  very  precise  tuning  of  their  opti¬ 
cal  path  and  therefore  are  susceptible  to  vibration  and  noise  amplification.  These 
disadvantages  make  FTSs  less  than  ideal  for  capturing  transient  scenes. 
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Figure  2.2.  Data  collected  from  FTS.  Results  are  then  converted  into  hypercubes  using 
Fourier  Transforms.  [12] 

2.3.2  Dispersive  Spectrometers. 

Diffraction  gratings,  prisms,  or  optical  filters  are  the  dispersive  components  often 
used  in  these  spectrometers.  Three  classes  of  dispersive  spectrometers  are  considered 
for  capturing  transient  events:  scanned-slit,  optical  filter,  and  tomographic.  Each 
class  collects  hypercubes  using  different  methods  schematically  represented  in  Fig¬ 
ure  2.3. 


Figure  2.3.  Example  of  the  spectrum  obtained  from  three  classes  of  dispersive  HSI 
devices.  Scanned-slit  and  filter  classes  require  scanning  while  the  tomographic  devices 
do  not.  [33] 
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The  x-y  plane  represents  the  FPA  and  the  red  band  is  the  spectral-spatial  infor¬ 
mation  collected  at  each  measurement.  Both  the  scanned  and  filter  classes  require 
multiple  collections  to  capture  the  entire  hypercube,  while  tomographic  instruments 
capture  the  hypercube  in  one  measurement.  Each  class  is  explained  in  the  next  sec¬ 
tions. 


2.3. 2.1  Scanned-Slit  Hyperspectral  Imagers. 

Scanned-slit  hyperspectral  imaging  systems  operate  by  dispersing  the  light  from 
a  slit  onto  the  FPA.  The  slit  is  narrowly-sized  to  achieve  high  spectral  resolution 
and  is  imaged  onto  one  dimension  of  the  FPA  which  is  the  x-axis  of  Figure  2.3.  The 
dispersed  light  is  spread  over  the  other  dimension,  the  y- axis,  of  the  FPA.  Figure  2.4 
shows  the  concept  of  a  conventional  scanned-slit  HSI. 


Figure  2.4.  Concept  of  a  scanned-slit  HSI.  The  system  images  a  slit  and  a  dispersive 
element  spreads  the  spectrum  along  one  dimension  of  the  FPA.  [12] 

To  obtain  a  hypercube,  these  devices  must  scan  for  the  remaining  spatial  dimen¬ 
sion.  These  instruments  use  pushbroom  or  wiskbroom  schemes  aboard  a  moving 
platform,  such  as  a  satellite  or  aircraft,  as  the  scanning  mechanism.  The  use  of  a 
slit  limits  the  area  coverage  at  any  given  time  and  reduces  the  optical  throughput 
therefore  reducing  the  signal-to- noise  ratio  (SNR).  Increasing  the  exposure  interval 
remediates  the  lack  of  signal  at  the  cost  of  system  response.  Both  of  these  limitations 
restrict  the  ability  to  capture  transient  events. 
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2.3. 2. 2  Spectral  Filters  Hyperspectral  Imagers. 


Spectral  filters  capture  two  spatial  dimensions  and  use  a  series  of  spectral  bandpass 
filters  to  obtain  the  spectral  content.  The  spectral  performance  is  limited  by  the 
bandpass  filters  coverage.  This  method  provides  good  area  coverage  and  the  ability 
to  spatially  capture  the  transient  event.  However,  it  must  sweep  through  the  spectral 
filters  therefore  cannot  capture  all  the  spectral  information  of  transient  events. 

2.4  Chromotomographic  Imagers 

CTIs  were  first  conceptualized  by  Russian  scientists  [13].  Since  its  concept,  a 
significant  amount  of  research  focusing  on  the  generation  of  proof-of-concept  devices 
[11,  32],  See  Mantravadi  and  Cain  [30]  for  a  more  extensive  history  of  chromotomo¬ 
graphic  imagers.  AFIT  began  working  on  modeling  chromotomographic  instruments 
in  2004  and  have  since  then  constructed  the  GCTEx  based  off  the  design  proposed 
by  Mooney  [31]. 

Unlike  conventional  HSIs,  which  require  scanning  or  filtering  in  either  the  spatial  or 
spectral  dimensions  to  generate  hypercubes,  chromotomographic  imagers  capture  all 
three  dimensions  (two  spatial  dimensions  and  one  spectral)  simultaneously.  Chromo- 
tomography  use  a  dispersive  element,  either  a  prism  or  diffraction  grating,  to  project 
the  three-dimensional  object  cube  onto  the  two-dimensional  FPA.  The  results  of  the 
projections  are  illustrated  in  Figure  2.5. 

As  shown  in  Figure  2.5,  the  image  formed  at  one  wavelength  is  multiplexed  with 
the  images  formed  at  all  the  wavelengths  over  the  spectral  range  of  the  dispersive 
element.  All  information  of  the  object  cube  is  captured  simultaneously,  no  loss  of 
information  due  to  scanning  or  filtering  occurs,  making  chromotomographic  imagers 
ideal  for  capturing  transient  events.  In  order  to  generate  hypercubes  of  the  object, 
projections  from  multiple  angles  are  required.  Furthermore,  a  complex  reconstruction 
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Figure  2.5.  Chromotomography  takes  a  3D  object  cube  and  projects  the  information 
onto  a  2D  detector  array.  The  varying  colors  represent  each  spectral  bin  of  the  object 
cube.  This  figure  illustrates  a  diffraction-based  chromotomgraphic  imager  where  a  set 
number  of  projections  are  spread  spectrally.  In  the  center  of  the  detector,  the  object 
cube  is  not  dispersed. 


algorithm  is  necessary.  Chromotomographic  imagers  offer  the  following  advantages 
over  conventional  HSIs: 


•  no  slit  or  filtering  required  therefore  high  throughput  enabling  lowered  SNR 
constraints 

•  high  temporal  resolution  due  to  the  lack  of  a  scanning  mechanism 

Diffraction  grating  chromotomographic  imagers  have  been  successfully  demon¬ 
strated  [11,  19].  They  offer  the  advantage  of  being  able  to  provide  multiple  pro¬ 
jections  instantaneously  which  may  be  advantageous  for  capturing  transient  events. 
However,  Mantravadi  demonstrated  in  simulated  scenarios  that  prism-based  chro¬ 
motomographic  imagers  offer  higher  spectral  resolution  than  diffraction-based  ones 
[29].  High  spectral  resolution  is  required  for  classification  of  detonation  events  [4], 
Thus,  the  focus  of  this  thesis  is  on  chromotomographic  imagers  that  use  prisms  as 
the  dispersive  element. 
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2.4.1  Principles  of  Operation. 


Chromotomographic  instruments  operate  under  similar  principles  as  medical  com¬ 
puted  tomography  (CT)  instruments.  Like  medical  CT  instruments,  projections  col¬ 
lected  at  multiple  angles  are  used  to  reconstruct  a  three-dimensional  cube  from  two- 
dimensional  detectors.  However,  the  third  dimension  in  chromotomography  cubes  is 
spectral  rather  than  spatial  [13]. 

Figure  2.6  shows  the  optical  layout  of  the  CTI  which  uses  a  rotating  DVP  to 
achieve  varying  spectral  dispersion.  As  shown  in  the  figure,  the  object  is  imaged 
through  the  optical  system  and  results  in  a  3-channel  hypercube.  Figure  2.6  demon¬ 
strates  that  the  hypercube  has  been  separated  into  its  red,  green,  and  blue  compo¬ 
nents. 


Figure  2.6.  Optical  layout  of  the  CTI  system  using  a  rotating  Direct  Vision  Prism  to 
vary  the  spectral  dispersion. 

The  telescope  segment  of  the  CTI  consists  of  two  lenses,  Li  and  L2,  spaced  to  the 
sum  of  each  focal  length  thereby  by  making  it  afocal.  Assuming  that  the  incoming 
light  into  the  first  lens  is  collimated,  the  light  refracted  by  the  second  lens  once  again 
becomes  collimated  and  angular  magnification  is  produced.  A  Field  Stop  (FS)  is 
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placed  on  the  intermediate  image  plane  at  the  focal  length  of  L3.  Since  the  prism 
spreads  the  object,  a  field  stop  is  required  to  ensure  the  imaged  object  does  not  exceed 
the  dimensions  of  the  FPA.  The  collimated  light  is  then  dispersed  by  the  prism  and 
then  is  imaged  by  the  focusing  lens,  L3,  onto  the  FPA. 

The  prism  rotates  the  spectral  contents  about  the  undeviated  wavelength,  shown 
in  green.  With  each  integration  time,  a  single  object  cube  projection  is  collected. 
Mooney  [31]  showed  that  the  radial  shift  from  the  center  of  the  FPA,  referred  to  as 
the  spectrometer  constant,  k,  is  defined  as 

k(X)  —  f3  tan  7(A),  (2.11) 

where  f3  is  the  focal  length  of  the  focusing  lens  and  7  is  the  dispersion  of  the  prism 
as  a  function  of  wavelength. 

Since  the  collected  data  has  spatial  and  spectral  information  interlaced  with  one 
another,  a  reconstruction  algorithm  is  required  to  deconvolve  the  spectral  information. 
The  following  section  discusses  the  reconstruction  algorithm. 

2.5  Reconstruction 

The  single  greatest  deficiency  in  chromotomographic  imager’s  performance  is  caused 
by  the  limited-angle  problem.  Due  to  the  nature  of  device  operation  coupled  with  the 
central  slice  theorem,  a  missing  cone  of  information  is  formed  preventing  all  of  the 
object  cube  from  being  collected.  In  order  to  produce  reliable  results,  reconstruction 
algorithms  must  fill  in  this  missing  information.  Therefore,  it  is  imperative  that  a 
full  understanding  of  the  central  slice  theorem  and  missing  cone  be  developed.  This 
section  is  then  followed  by  the  analytical  expressions  that  are  common  amongst  the 
reconstruction  algorithms . 
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2.5.1  Central  Slice  Theorem. 


The  Fourier  slice  theorem,  commonly  referred  to  as  the  central  slice  theorem,  is 
the  fundamental  theory  that  governs  tomographic  systems.  The  central  slice  theorem 
is  defined  as  [24]: 

Definition.  The  Fourier  transform  of  a  parallel  projection  of  an  object  f(x,  y)  ob¬ 
tained  at  angle  6  equals  a  line  in  a  2D  Fourier  transform  of  f(x,  y)  taken  at  the  same 
angle. 

By  performing  the  Fourier  transform  on  a  projection,  a  line  through  the  center  in 
the  2D  Fourier  transform  of  the  object  is  obtained.  Figure  2.7  is  an  example  of  the 
central  slice  theorem  of  a  image  collected  by  a  medical  CT  instrument. 


p(x)  =  £  f(x,y)dy  P(u)=  £  £  f{x,y)e  l2mixdxdy  P(u)=  £  £  f{x,y)e  l2,ruxdxdy 

J—  oo  d—  GO  d—  GO  00  J—  00 


Figure  2.7.  Illustration  of  central  slice  theorem  for  a  single  projection  collected.  The 
Fourier  transform  of  the  projection  is  equal  to  a  line  through  the  center  in  the  2D 
Fourier  transform  of  the  object.  [24] 


If  a  sufficient  number  of  projections  is  collected  over  the  range  from  0  to  n,  the 
entire  Fourier  space  of  the  object  being  reconstructed  is  filled.  Once  the  object’s 
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Fourier  transform  is  obtained,  the  object  can  be  fully  recovered  by  taking  the  inverse 
Fourier  transform.  The  tomographic  reconstruction  process  is  a  series  of  ID  Fourier 
transforms  followed  by  a  2D  inverse  Fourier  transform. 

2.5.2  The  Missing  Cone. 

The  central  slice  theorem  extends  to  three  dimensions.  For  a  object  cube,  the  2D 
Fourier  transform  of  the  tomographic  projection  is  equal  to  the  3D  Fourier  transform 
of  the  image  evaluated  on  a  plane  through  the  origin  in  a  direction  perpendicular  to 
the  projection  bean  [6],  this  is  called  the  x-ray  projection  plane.  The  projection  beam 
angle  is  controlled  by  the  spectrometer  constant  defined  in  Equation  (2.11).  Figure  2.8 
depicts  how  the  object  cube  creates  tomographic  projections  onto  the  detector  via 
the  x-ray  transform.  The  recovery  of  a  3D  distribution  from  2D  projections  is  known 
as  the  x-ray  transform  [11], 

Similar  to  the  2D  CT  reconstruction  discussed  in  Section  2.5.1,  the  entire  3D 
Fourier  space  of  the  object  needs  to  be  filled  to  reconstruct  the  object  cube.  By 
rotating  the  prism,  different  slices  through  the  origin  are  obtained.  However,  the 
same  projection  beam  angle  remains  constant  throughout  all  the  prism  rotations. 
The  constant  projection  beam  and  failure  to  fill  in  the  entire  3D  Fourier  space  is 
known  as  the  limited-angle  problem  for  tomographic  instruments. 

In  the  limit  of  an  infinite  number  of  different  projections,  the  projections  will  define 
a  cone,  where  the  half-angle  of  the  cone  vertex  is  approximated  to  be  45  degrees  [32], 
Since  none  of  the  projections  provide  any  information  about  the  image  inside  this 
cone,  it  is  referred  to  as  the  missing  cone.  Figure  2.9  illustrates  how  the  missing 
cone’s  shape  is  formed  in  the  frequency  space. 

As  clearly  evident,  there  are  gaps  within  the  frequency  space.  Discussed  in  the 
next  section  in  more  detail,  these  gaps  in  the  frequency  space  are  the  reason  why 
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Figure  2.8.  Geometry  of  chromotomographic  data  collection  and  its  relation  to  the 
x-ray  transform.  [6] 


a  direct  inversion  method  for  producing  a  hypercube  is  ill-posed.  Some  mechani¬ 
cal  techniques  of  filling  in  the  missing  cone  described  in  [32]  exist;  however,  some 
reconstruction  algorithms  is  devised  to  fill  in  the  missing  cone.  It  is  worth  noting 
that  a  monochromatic  source  would  fill  in  all  the  frequency  space  information  as  a 
monochromatic  object  is  essentially  a  two-dimensional  object,  reverting  back  to  the 
2D  central  slice  theorem. 

2.5.3  System  Transfer  Function. 

The  System  Transfer  Function  (STF)  is  shared  among  all  the  reconstructions  in 
literature.  The  STF  is  a  matrix  that  contains  the  location  of  all  the  spectral  point 
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Figure  2.9.  2D  Fourier  transform  of  each  projection  equals  a  plane  through  the  3D 
frequency-space^ ,  £2»  C)  representation  of  the  object  cube.  The  missing  cone  is  the 
result  of  limited- angle  problem  and  the  inability  to  fill  in  all  of  the  frequency  space.  As 
the  number  of  projections  increases,  as  shown  in  (a)-(c),  the  shape  of  the  missing  cone 
becomes  more  prominent.  [14] 


spread  functions  created  from  the  prism. 

Mooney  showed  in  [31]  that  the  measured  projected  data,  g(xi,x2,(f>),  is  found 
by  summing  each  tomographic  projection  of  the  desired  hypercube,  f(x ±,x2,  A)  with 
respect  to  A 

/OO 

f(x i,  x2,  A)  *  —  fc(A)  cos(0) ,  x2  -  k( A)  sin(0))dA,  (2.12) 

-OO 


where  X\  and  x2  are  the  pixel  coordinates  of  the  measured  data  cube  and  hypercube, 
0  is  the  projection  angle  bound  between  0  and  27 r,  and  k( A)  is  the  spectrometer 
constant  defined  in  Equation  (2.11).  Similar  to  the  imaging  convolution  integral 
in  Equation  (2.1),  the  imaged  formed  is  the  result  of  shifts  from  prism  onto  the 
object.  These  shifts  are  summed  for  all  wavelengths  yielding  the  measured  data 
cube.  Equation  (2.12)  can  be  expressed  more  succinctly  as 

/OO 

f(x-  k(X)p<j),X)dX,  (2.13) 

-OO 

where  x  =  (xi,x2),  and  p $  =  (cos(0),  sin(0)).  Taking  the  2D  Fourier  transform  of 
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Equation  (2.13)  becomes 


G(£  <f>) 


exp[-2rt(A)£  ■  A)dA, 


(2.14) 


where  £  is  the  spatial  frequency  representation  of  x ,  and  £-p^  denotes  the  vector  dot 
product.  Considering  that  the  data  cubes  are  sampled  at  a  discrete  spectral  band,  n, 
and  discrete  projection,  m,  Equation  (2.14)  is  now  a  summation 


N- 1 

Gm(0  =  ^2  exP[-27rn£  •  Pm]K(€,  A)dA,  (2.15) 

n=0 

where  pm  =  (cos(2-7rm/M), sin(27rm/M))  for  M  projections  measured,  0  <m<  M, 
and  n  =  k( A)  for  spectral  bands  of  the  device. 

In  terms  of  linear  algebra,  the  spatial  frequency  representations  of  the  projection 
images  and  object  cube  are  related  by 

G(0  =  A(0F(0,  (2.16) 

where  A(£)  is  a  M  x  N  matrix  known  as  the  STF.  For  each  spatial  frequency,  the 
elements  of  A(£)  are 

An,n(£)  =  exp[-27rn£  •  pm\.  (2.17) 

The  image  cube  is  obtained  by  inverting  7l(£), 

F(0  =  AtiV'G®.  (2.18) 

In  order  to  invert  the  matrix,  M  >  N ;  however,  A(£)  is  often  rectangular  and  does 
not  have  full  rank.  Singular  Value  Decomposition  (SVD)  is  used  to  find  the  pseudo¬ 
inverse  and  the  minimum-norm  solution  is  obtained.  Since  there  are  gaps  within  the 
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frequency  space,  advanced  reconstruction  algorithms  must  fill  in  the  cone  of  missing 
information  to  achieve  more  accurate  results. 

2.6  Reconstruction  Development 

An  alternative  approach  in  finding  the  minimum-norm  solution  is  to  use  a  filtered 
backprojection  (FBP)  technique  modeled  from  the  medical  tomography  field  [2],  Fil¬ 
tered  backprojections  essentially  “shift- and- add”  the  collected  projection  images.  For 
each  projection  angle,  m,  of  the  spectral  bin,  n,  the  projected  image  is  shifted  by  the 
conjugate  of  Equation  (2.17).  This  centers  the  projected  imaged  for  the  given  spectral 
bin.  For  all  M  projections,  the  shifted  projection  images  are  summed  resulting  in  a 
“stack”  of  shifted  projected  images.  As  the  shifted  projected  images  are  summed, 
the  information  within  the  spectral  bin  increases  its  structural  integrity  and  over¬ 
come  the  residual  artifacts  from  the  other  spectral  bin  content.  Filtering  mechanisms 
are  applied  at  either  each  successive  shift  or  after  summing  the  shifted  images  to  re¬ 
duce  the  residual  artifacts.  The  filtered  backprojection  algorithm  needs  one  iteration 
and  is  computationally  efficient  as  Fourier  transforms  and  matrix  multiplications  are 
required. 

To  fill  in  the  missing  cone,  an  iterative  process  such  as  Projection  Onto  Convex 
Sets  (POCS)  or  maximum  likelihood  algorithms  is  required.  The  iterative  recon¬ 
struction  algorithms  achieve  recovery  of  the  unknown  object  by  utilizing  a  priori 
information  about  the  image.  Typically,  an  initial  estimate  about  the  unknown  ob¬ 
ject  is  made,  which  is  then  subjected  to  a  sequence  of  corrections,  forcing  it  to  satisfy 
a  number  of  desirable  characteristics.  This  sequence  of  corrections  is  applied  repeti¬ 
tively  until  convergence  occurs.  Mooney  et  al.  developed  algorithms  that  used  SVD 
to  calculate  the  minimum-norm  solution  and  feed  it  into  a  iterative  POCS  [6,  7]  algo¬ 
rithm.  Gould  developed  an  reconstruction  algorithm  using  Estimation  Theory  that 
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uses  the  STF  to  produce  the  maximum  likelihood  estimate  of  the  object  [16].  Both 
of  these  algorithms  demonstrated  the  ability  to  reconstruct  the  scene  accurately. 

A  computationally  efficient  reconstruction  algorithm  is  desirable  for  AFIT.  The 
space-based  CTI  has  embedded  processing  and  runs  the  risk  of  overheating  [27]. 
Lowering  the  computational  demands  of  the  algorithm  reduces  the  risk  of  overheating. 
In  that  spirit,  a  computationally  efficient  filtered  backprojection  algorithm  developed 
by  Deming  in  [10]  is  implemented  for  this  thesis. 

2.6.1  Fast  Reconstruction  Algorithm. 

The  fast  reconstruction  algorithm  is  a  filtered  backprojection  which  obtains  the 
minimum-norm  solution  using  closed-form  expressions.  The  algorithm  avoids  calcu¬ 
lating  the  pseudo- inverse  by  using  a  series  of  matrix  multiplications,  Fast  Fourier 
Transforms  (FFTs),  and  applying  a  spectral  filter.  Fewer  computations  are  needed 
in  calculating  the  FFT  verse  computing  the  pseudo-inverse.  The  amount  of  process¬ 
ing  required  is  compounded  when  dealing  with  large  focal  plane  array  formats  and 
collecting  a  large  number  of  projections. 

The  reconstruction  algorithm  can  be  summarized  as  follows.  For  each  spatial 
frequency,  the  preliminary  (unfiltered)  Fourier  domain  representation  of  the  image 
cube  (referred  to  as  the  Fourier  cube),  F'(£,  A),  is  computed  by  first  backprojecting 
the  data,  G(£,  4>)  with  the  Hermitian  adjoint  of  A(£) 

/*2tt 

F'(ZA)  =  [AHG](Z,\)=  exp[2nik(\)£-p,i,\G(Z,<i))d<f>.  (2.19) 

Jo 

The  filtered  Fourier  cube,  F/x(^,  <U\),  is  found  by  ID  Fourier  transforming  the 
unfiltered  Fourier  cube  from  Equation  (2.19)  relative  to  A  to  compute  F'(£,u> \)  and 
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then  applying  the  spectral  filter 


F,(^ux) 


a) 

ji  +  M(£,u;  a) 


(2.20) 


The  filtering  parameter,  n,  and  the  spectral  filter,  M(£,uj\),  are  derived  from  filtering 
methods  in  the  medical  CT  FBPs.  M (£,  u\)  is  a  zero-order  Bessel  function  that  has 
been  ID  Fourier  transformed  with  respect  to  A,  simplified  to 


M(t,ux)=2(k2er-uiri/2 


(2.21) 


where  =  (£1 +£f)1/2  and  k  is  the  spectrometer  constant  defined  in  Equation  (2.11). 
The  final  step  is  inverse  Fourier  transforming  with  respect  to  c 0\  to  obtain  the  filtered 
Fourier  cube.  The  value  of  fi  depends  on  the  noise  level  in  the  data,  ideally  ft  will 
approach  0  but  remain  large  enough  to  ensure  stability  in  the  presence  of  noise. 


2.7  Spatial  and  Spectral  Resolution  Limits 

The  limits  of  the  spatial  and  spectral  resolution  is  based  upon  the  most  sim¬ 
plistic  of  objects,  a  point-source.  In  [5],  Bostick  built  mathematical  models  and 
collected  experimental  data  using  point-sources  to  determine  the  spatial  and  spec¬ 
tral  resolution.  Bostick  concluded  that  the  spatial  resolution  is  defined  by  full- width 
half-maximum  (FWHM)  of  the  point  spread  function  as  asserted  by  the  Rayleigh  Cri¬ 
terion.  If  the  optical  system  is  aberration  free,  the  spatial  resolution  will  be  diffraction 
limited.  The  spatial  features  should  not  be  impacted  as  the  result  of  the  prism  disper¬ 
sion,  this  is  evident  when  considering  a  monochromatic  object.  The  spatial  features 
of  a  monochromatic  object  are  preserved  as  there  is  no  other  spectral  information 
is  multiplexed  with  it.  In  the  presence  of  aberrations,  the  size  of  the  PSF  will  blur 
thereby  limiting  the  ability  to  resolve. 


Bostick  determined  that  the  limit  of  the  spectral  resolution  is  dependent  on  the 
prism  dispersion  and  the  FWHM  of  the  spectral  PSF.  Under  the  assumption  the 
reconstruction  algorithm  works  ideally,  if  two  objects  which  require  a  minimum  spatial 
separation  of  w  to  be  resolved,  the  minimum  spectral  separation  AAm  given  by 
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(2.22) 


In  other  words,  the  change  in  radial  shift  from  one  wavelength  to  another  must  be 
greater  than  the  FWHM  of  the  spot  size  for  the  wavelengths  of  interest.  Consider 
a  polychromatic  point  source  in  a  diffraction-limited  CTI  system  and  ignoring  per¬ 
formance  and  operation  of  the  FPA.  In  order  to  resolve  the  spectral  contents,  the 
radial  shift  from  the  prism  must  be  greater  than  the  radius  of  the  Airy  disk  given  by 
Equation  (2.8).  If  aberrations  are  present,  the  size  of  the  Airy  disk  increases  thereby 
lowering  the  spectral  resolution. 


2.8  Ground-based  CTEx 

AFIT  has  developed  several  iterations  of  laboratory  and  held  instruments,  col¬ 
lectively  referred  to  as  GCTEx,  over  the  past  several  years.  With  each  revision,  the 
devices  performance  and  capabilities  have  incrementally  increased.  The  first  lab  tests 
where  conducted  by  LeMaster  [28]  followed  by  Bostick  [5].  O’Dell  [36]  created  a 
Newtonian  based  instrument  and  collected  the  Erst  set  of  held  data.  This  system  was 
challenging  to  align  and  did  not  maintain  calibration  which  made  deployment  into  the 
held  difficult.  Niederhauser  [35]  revamped  the  GCTEx  into  a  linear  system  where  all 
the  components  are  centered  along  the  optical  axis,  resulting  in  better  performance 
and  more  robust  held  deployment. 
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2.8.1  Optical  System. 


The  linear  system  developed  by  Niederhauser  [35]  follows  the  layout  of  Figure  2.6. 
It  utilizes  three  telephoto  lenses  to  mitigate  aberrations,  specified  in  Table  2.4,  and 
the  spinning  DVP  is  housed  in  a  solid  block  of  aluminum  to  reduce  vibrations  from 
the  motor.  The  linear  system  is  machined  such  that  no  adjustments  are  necessary 
and  is  quick  to  assemble. 


Table  2.4.  GCTEx  Lens  Specifications. 


Lens 

Model 

Focal  Length  (mm) 

F/# 

Mount  Type 

U 

AF-S  Nikkor  400m  f/2.8G  ED  VR 

400 

2.8 

F 

L2 

Tarnron  1A1HB  75mm 

75 

3.9 

C 

L3 

AF  Nikkor  85mm  f/1.8D 

85 

1.8 

F 

The  linear  system  is  shown  in  Figure  2.10.  The  system  uses  a  high-speed  video 
camera,  Phantom  v5.1,  from  Vision  Research  which  collects  images  at  rates  up  to  10 
kHz. 


Figure  2.10.  Linear  GCTEx  instrument  developed  by  Niederhauser.  [35] 


The  DVP  is  housed  in  a  circular  shaft  that  spins  at  rate  up  to  25  rev/s.  For  each 
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revolution,  a  complete  hypercube  produced,  therefore  the  GCTEx  produces  up  to  25 
hypercubes  per  second.  For  the  space-based  system,  the  prism  will  rotate  nominally 
at  100  rev/s.  The  DVP  was  designed  by  Bostick  and  is  made  of  an  optically-bonded 
assembly  of  Schott  LaSF  N30  and  SFL6  glass  8  components  as  shown  in  Figure  2.11. 


Figure  2.11.  Geometry  and  components  of  the  Direct  Vision  Prism.  [36] 


The  dispersion  of  the  prism  could  be  solved  using  Snell’s  law;  however,  it  has  been 
more  convenient  to  express  the  deviation  angle  caused  by  the  prism  in  the  form  of 
the  power  model, 

7  =  a\b  +  c.  (2.23) 

The  DVP  curve  fit  parameters:  a,  b,  c  for  the  deviation  angle  (degrees),  7,  for  a  given 
wavelength  (/tm)  are  listed  in  Table  2.5.  The  data  used  for  the  power  model  is  based 
off  of  calculations  from  Zernax  model  of  the  DVP.  The  undeviated  wavelength,  7  =  0, 
of  the  prism  is  at  ^547  nm. 

The  linear  system  demonstrated  that  the  optical  performance  increased  from  the 
Newtonian  design.  The  polychromatic  MTFs  measured  where  significantly  higher 
and  most  importantly,  the  measured  deviation  angles  matched  the  theoretical  values 
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Table  2.5.  Theoretical  DVP  Curve-fit  Parameters  from  by  Zemax  model. 

Term  Value 


a  .2182 

b  -3.329 

c  -1.633 


predicted  by  Zemax.  In  Figure  2.12,  three  emission  lines  from  a  mercury  pen  lamp 
match  the  Zemax  model. 
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Figure  2.12.  Mercury  pen  lamp  emission  lines  matching  the  Zemax  model.  [35] 


The  match  in  deviation  is  crucial  for  reconstruction;  however,  tests  of  the  linear 
system  concluded  that  a  transverse  offset  which  laterally  shifts  the  deviation  angle 
from  its  predicted  center  was  present.  Additionally,  since  the  GCTEx  uses  telephoto 
lenses  where  multiple  groups  of  lenses  are  used  to  reduce  monochromatic  aberrations, 
the  system  has  a  significant  amount  of  chromatic  aberrations.  Both  of  chromatic 
aberrations  and  transverse  offset  are  discussed  in  the  next  sections. 
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2.8. 1.1  Chromatic  Aberrations. 


To  demonstrate  chromatic  aberrations,  mercury  (Hg)  and  neon  (Ne)  pen  lamps 
are  used  in  the  laboratory.  Mercury  and  neon  pen  lamps  are  used  to  calibration 
sources  because  of  of  their  distinct  spectral  lines  features.  Table  2.6  lists  the  spectral 
lines  and  the  corresponding  relative  normalized  intensity  [34,  38]  for  each  source  used 
throughout  this  thesis. 


Table  2.6.  Normalized  Relative  Intensity  Emission  Lines  of  Hg  and  Ne  Pen  Lamps. 


Hg  Lines  /im 

Ne  Lines  /.irn 

Normalized  Relative  Intensity 

0.40465 

0.44 

0.43583 

1.00 

0.54683 

1.00 

0.57696 

0.11 

0.57907 

0.12 

0.5853 

0.2000 

0.5882 

0.1000 

0.5945 

0.0500 

0.6096 

0.0300 

0.6153 

0.1000 

0.6217 

0.1000 

0.6267 

0.1000 

0.6334 

0.1000 

0.6402 

0.2000 

0.6506 

0.1500 

0.6601 

0.0100 

0.6678 

0.0500 

0.6717 

0.0070 

0.6923 

1.0000 

0.7028 

0.3400 

0.7245 

0.7700 

0.7489 

0.3200 

A  point  source  was  simulated  by  placing  and  iris  in  front  of  the  mercury  and 
neon  pen  lamps.  The  results  taken  using  GCTEx  are  shown  in  Figures  2.13(a)  and 
Figure  2.13(b). 
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(a)  (b) 


Figure  2.13.  (a)  Chromatic  blurring  of  the  435.8  nm  line,  near  the  top  of  the  image, 
from  a  Hg  pen  lamp,  (b)  The  shorter  wavelengths,  closest  to  the  top  of  the  image,  of 
the  Ne  pen  lamp  are  in  focus  while  the  longer  wavelengths  are  blurred. 

From  the  mercury  pen  lamp,  the  FWHM  of  the  546.8  nm  line  is  2.30  pixels  while 
the  435.8  nm  line  has  a  FWHM  of  7.37  pixels. 

2. 8. 1.2  Transverse  Offset. 

The  transverse  offset,  or  pinwheel  offset,  presents  a  serious  problem  for  the  recon¬ 
struction  algorithm.  The  fast  reconstruction  algorithm  relies  on  the  location  of  the 
spectral  PSFs.  Not  accounting  for  the  offset  results  in  an  error  in  the  location  and 
thus  a  significant  performance  drop  in  reconstructing  the  scene.  Figure  2.14  shows 
measured  center  is  offset  from  the  predicted  center  of  the  circular  dispersion  path. 
It  is  important  to  note  that  dispersion  within  the  radial  line  remains  consistent  with 
the  predicted  dispersion  values. 

Niederhauser  equated  the  offset  as  a  jump  discontinuity  and  averaged  the  radii 
to  obtain  a  correction  in  the  power  model  fit  parameters.  However,  it  is  important 
to  realize  the  lateral  shift  is  constant  for  all  wavelengths  and  is  symmetric  for  all 
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Figure  2.14.  Demonstration  of  the  transverse  offset  from  a  mercury  pen  lamp  point 
source.  The  dispersion  along  the  “true”  deviation  line  matches  the  predicted  model; 
however,  the  transverse,  or  “pinwheel’,  offset  shifts  the  center  of  the  deviation  from 
the  predicted  center.  [35] 


rotations.  Figure  2.15  displays  a  mercury  pen  lamp  point  source  at  projection  angles 
in  90  degree  increments.  For  each  projection  angle,  the  435.8  nm,  546.8  nm,  and 
576.9  nm  emission  lines  are  shown  with  the  435. 8nm  line  being  the  outermost. 

For  the  546.8  nm  line,  a  blue  circle  is  added  to  illustrate  that  offset  is  radially  sym¬ 
metric  over  the  projections.  Also  note  there  is  a  phase  lag  for  each  projection.  Ideally, 
all  of  the  projections  for  the  546.8  nm  line  would  align  perfectly  with  dashed  green 
line  which  is  the  center  of  the  radial  offset.  The  phase  lags  for  shorter  wavelengths 
under  the  undeviated  wavelength  and  leads  for  the  longer  wavelengths.  This  phase 
lag  changes  with  wavelength  and  must  also  be  accounted  for  in  the  reconstruction 
algorithm. 
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Figure  2.15.  Four  orthogonal  projections  measured  using  a  mercury  pen  lamp  point 
source. 

2.8.2  Phantom  Camera. 

GCTEx  uses  the  Vision  Research  Phantom  v5.1  monochrome  CMOS  sensor.  The 
camera  is  capable  of  collecting  data  at  fast  frame  rates,  up  to  95  kHz  at  the  smallest 
resolution  of  64x32  pixels.  However,  the  dispersion  of  the  prism  and  a  85  mm  focus¬ 
ing  lens  require  that  512x512  resolution  be  used,  yielding  a  maximum  frame  rate  of 
4380  frames  per  second.  If  larger  formats  are  required,  the  maximum  rate  will  drop 
accordingly.  Table  2.7  summarizes  the  list  of  camera  specifications.  The  pixels  of  the 
camera  are  assumed  to  be  square. 

The  camera  is  able  to  frame  quickly  by  storing  the  collected  data  to  internal 
memory.  The  camera  will  continuously  run  until  a  trigger  is  received  by  the  user. 
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Table  2.7.  Phanton  v5.1  Camera  Spefications. 


Parameter 

Value 

Detector 

Silicon 

Full  Resolution  (pixels) 

1024x1024 

Frame  Rate  at  Full  Resolution  (Hz) 

1200 

Minimum  Exposure  Time  (/xs) 

2 

Pixel  Pitch  (gm) 

20 

Bit  Depth  (bits) 

8 

On-board  Memory(MBytes) 

4096 

The  frames  are  placed  in  a  circular  buffer,  so  the  oldest  frame  is  overwritten  once 
the  memory  is  full.  The  camera’s  system  gain,  noise  performance,  and  quantum  effi¬ 
ciency  (QE)  are  necessary  for  estimating  the  amount  of  signal  necessary  to  detect  a 
target.  The  manufacturer  provides  documentation  on  the  camera  parameters;  how¬ 
ever,  the  documentation  is  limited  and  not  sufficient.  A  photon  transfer  curve,  as 
described  in  [25],  was  generated  to  determine  system  gain  and  noise  performance. 
A  monochromator  and  calibrated  photodiode  setup  were  used  to  calculate  quantum 
efficiency.  The  results  from  the  photon  transfer  curve  are  listed  in  Table  2.8. 


Table  2.8.  System  gain,  KadCi  and  read  noise  of  the  Phantom  Camera. 


Parameter 

Value 

Uncertainty 

Kadc  (e"/DN) 

143.12 

±  7.73 

Read  noise  (e_) 

34.87 

±  7.92 

The  camera’s  performance  is  poor,  the  sensitivity,  Kadc ,  of  the  camera  is  low 
and  the  amount  of  photo-electrons  are  required  to  rise  above  read  noise  floor  is  high. 
Furthermore,  the  camera  has  considerable  amounts  of  fixed-pattern  noise  (FPN), 
the  pixcl-to-pixel  gain  non-uniformity,  which  had  to  be  subtracted  to  obtain  these 
results.  Additionally,  the  camera  exhibited  a  nonlinearity  in  the  system  gain.  This 
nonlinearity  prevents  a  flat-field  correction  to  be  applied  over  the  entire  dynamic 
range.  The  poor  performance  and  the  nonlinearity  are  thought  to  be  attributed  to 
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the  high  speed  nature  of  the  device.  These  results  are  reported  but  ignored  and 
assumed  to  be  ideal  for  the  purposes  of  creating  a  simulation  and  evaluating  the 
reconstruction  algorithm. 

The  QE  of  the  camera  was  calculated  from  the  wavelength  from  400  nm  to  740 
nm  in  increments  of  20  nm.  The  QE  measurements  were  taken  at  lower  frame  rate 
in  order  to  increase  the  integration  time  to  achieve  a  higher  signal.  The  results  are 
shown  in  Figure  2.16  below. 
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Figure  2.16.  Quantum  efficiency  over  the  visible  regime.  The  oscillations  are  consistent 
with  those  observed  by  O’Dell  and  the  camera  manufacturer’s  datasheet. 

The  camera’s  QE  has  a  peak  response  at  460  nm  of  0.42  and  remains  above  0.20 
over  the  entire  wavelength  range.  As  expected,  the  shorter  wavelengths  have  lower 
QE  due  to  reflections.  At  the  longer  wavelengths,  QE  drops  off  because  it  is  more 
difficult  to  generate  photo-electrons  as  the  energies  approach  the  required  band  gap 
energy.  The  QE  measurements  published  by  the  camera  vendor  advertise  QE’s  above 
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0.50  from  450  nm  to  650  nm;  however,  these  measurements  published  were  for  a 
newer  model  [41].  The  QE  oscillates  over  the  spectral  range  and  the  oscillation  is 
more  pronounced  if  smaller  increments  are  used  and  is  averaged  out  with  20  nm 
increments.  The  shape  of  the  QE  curve  matches  the  relative  measurements  made  by 
O’Dell  [36]  and  are  consistent  with  datasheets  from  the  manufacturer. 

2.8.3  Electrical  Support  System. 

Fixing  the  electronics  system  was  one  of  the  main  research  goals  in  this  thesis. 
Knowledge  of  the  prism  angle  is  imperative  for  the  reconstruction  algorithm.  Previ¬ 
ous  AFIT  had  to  manually  calculate  prism  angles  based  off  of  poorly  synchronized 
timestamps  and  made  it  difficult  for  reconstructing  scenes.  Thus,  the  entire  elec¬ 
tronics  and  software  system  were  redeveloped  in  order  to  fix  the  logging  of  the  prism 
angles  in  relation  to  the  camera  files. 

The  GCTEx  Electrical  System  controls  the  speed  of  the  prism,  operation  of  the 
Phantom  v5.1  camera,  and  collects  the  encoder  data  associated  with  the  prism’s  angle. 
A  rackmount  computer  houses  the  SAM-3  data  acquisition  module  and  controls  the 
camera.  The  motor  control  and  a  custom  electronics  board  that  converts  the  encoder 
signals  into  prism  angles  are  housed  in  a  custom  enclosure.  Finally,  a  network  switch 
connects  the  camera,  computer,  and  any  other  peripheral  using  the  local  network. 
All  of  the  electrical  equipment  is  bundled  into  a  small  8U  portable  rack  for  ease  of 
transportation  as  shown  in  Figure  2.17.  See  Appendix  B  for  the  electrical  system 
block  diagram. 

The  prism  angle  is  measured  by  a  quadrature  encoder  and  has  a  resolution  of 
4096  bits  per  revolution.  This  means  that  the  angular  resolution  of  the  prism  angle  is 
within  0.088  degrees  or  1.534  milliradians.  The  prism  angle  is  sampled  by  the  SAM-3 
module  and  is  synchronized  with  the  captured  images  using  timing  signals  provided 
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Figure  2.17.  Portable  CTEx  rack  for  easy  transportation  to  field. 


by  the  camera.  The  timing  signals  allow  the  prism  angle  to  be  captured  at  the  start 
and  end  of  the  camera’s  integration  time.  This  architecture  avoids  interpolating  the 
prism  angles  using  timestamps  from  the  SAM-3  and  camera,  which  posed  a  significant 
problems  for  previous  AFIT  students. 

The  electronics  are  capable  of  sampling  the  encoder  data  at  1  MHz,  which  corre¬ 
lates  to  the  maximum  speed  of  244  rev/sec  before  the  prism  has  traveled  a  distance 
more  than  the  encoder  resolution.  This  maximum  speed  is  well  above  the  notional 
operating  speed  of  100  rev/sec. 
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3.  System  Model  Methodology  and  Validation 


The  need  for  a  high  fidelity  simulation  of  the  optical  system  is  imperative  for 
future  development  of  AFIT’s  chromotomography  experiment  (CTEx)  project  as  the 
system  has  yet  to  be  fully  characterized  and  understood.  Furthermore,  a  high  fidelity 
simulation  provides  an  invaluable  tool  for  developing  the  reconstruction  algorithm 
and  determining  design  criteria  for  the  spaced-based  system,  referred  to  as  SCTEx. 

CTEx’s  optical  system  has  been  modeled  in  Zemax  over  the  past  several  years. 
Zemax  is  an  industry  standard  for  computer-aided  optical  system  design  and  is  a 
well-vetted  optical  system  analysis  tool.  Zemax  is  a  physics-based  model  that  uses 
geometric  ray  traces  and  wave  propagation  to  obtain  more  realistic  simulations.  Ze- 
max’s  capabilities  will  be  heavily  leveraged  upon  to  build  the  simulation.  The  sim¬ 
ulation  model  uses  MATLAB  to  serve  as  an  interface  to  Zemax  using  the  Dynamic 
Data  Exchange  (DDE)  interface  tools  created  by  Griffith[17].  The  design  goals  for 
the  simulation  are: 

•  Provide  a  mechanism  to  simulate  a  variety  of  scenarios 

•  Move  away  from  paraxial  approximations  and  develop  more  representative  mod¬ 
els  that  demonstrate  the  effects  of  aberrations 

•  Decouple  the  simulation  from  the  model,  i.e.  the  simulation  will  work  with  any 
Zemax  model 

3.1  Assumptions 

Several  assumptions  are  made  to  reduce  complexity  while  maintaining  the  integrity 
of  the  simulation: 

1.  Transmission  loss  due  to  atmospheric  absorption  is  negligible  operating  in  the 
VIS/NIR 
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2.  No  loss  of  image  quality  or  transmission  dne  to  atmospheric  effects 

3.  All  sources  are  incoherent 

3.2  Overview  of  Simulation 

The  foundation  of  the  chromotomographic  simulation  lies  within  the  Zemax  model 
for  the  optical  system.  In  order  for  the  simulation  to  be  accurate,  requires  that 
the  Zemax  model  be  accurate  as  well.  Once  the  Zemax  model  is  developed,  the 
simulation  can  be  executed.  The  simulation  is  setup  by  MATLAB  so  interchanging 
Zemax  models  requires  minimal  effort.  This  capability  allows  different  Zemax  models 
to  be  quickly  evaluated. 

The  DDE  interface  allows  MATLAB  to  control  Zemax  functions  autonomously. 
Almost  all  Zemax  functions  are  available  through  the  DDE  interface,  if  they  are  not 
available,  functions  can  be  created  using  Zemax’s  scripting  language.  MATLAB  is 
able  to  rotate  the  prism,  vary  wavelengths,  and  perform  image  analyses.  Figure  3.1 
is  a  flowchart  of  the  simulation  using  all  the  DDE  function  calls. 

The  simulation  defines  the  object,  both  spatially  and  spectrally,  as  a  function 
of  prism  angle.  Time  is  controlled  by  varying  the  prism  speed,  and  the  angular 
displacement  of  the  prism  angles  varies  with  according  to  prism  speed.  For  each 
simulation,  the  user  defines  the  spatial  and  wavelength  profiles  of  interest. 

Using  the  set  spatial  profile,  MATLAB  creates  or  defines  an  image  source  hie 
to  be  simulated  at  each  prism  rotation.  An  image  simulation  is  completed  for  all 
wavelengths  listed  in  the  wavelength  profile.  The  results  of  the  image  simulation 
are  then  scaled  according  the  the  wavelength  profile.  After  all  wavelengths  have 
been  simulated  and  scaled,  the  results  of  all  simulations  are  summed  to  produce  the 
simulated  output  for  each  prism  angle  defined. 
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Figure  3.1.  Flowchart  of  the  chromotomographic  simulation. 
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3.2.1  Zemax  Image  Simulation. 


The  Zemax  image  simulation  is  based  upon  Equation  (2.2),  where  a  grid  of  point 
spread  functions  are  convolved  with  the  image  source  file  to  form  the  simulated  image 
through  the  optical  system.  The  grid  spans  the  field  size  and  describes  the  aberrations 
at  selected  points  in  the  field  of  view  defined  by  the  source  image  file  and  field  size 
settings. 

The  Image  Simulation  tool  performs  the  following  operations  to  produce  the  sim¬ 
ulated  image: 

1.  The  source  image  file  is  oversampled  and  a  guard  band  is  applied  specified  by 
the  user 

2.  A  grid  of  PSFs  are  computed 

3.  The  PSF  grid  is  interpolated  for  every  pixel  in  the  modified  source  image  file 

4.  At  each  pixel,  the  effective  PSF  is  convolved  with  the  modified  source  image 
file  to  determine  the  aberrated  bitmap  image 

5.  The  resulting  image  bitmap  is  then  scaled  and  stretched  to  account  for  the 
detected  image  pixel  size,  geometric  distortion,  and  lateral  color  aberrations 

An  example  of  the  image  simulation  PSF  grid  and  the  resulting  simulated  image 
are  shown  in  Figures  3.2(a)  and  3.2(b),  respectively. 

The  calculation  of  the  PSF  grid  is  critical  for  accurate  results.  There  are  two 
methods  for  modeling  the  aberrations  using  the  PSF:  diffraction  and  geometric.  The 
diffraction  method  computes  the  PSF  using  direct  integration  of  Huygens  wavelets 
method[44].  The  geometric  method  computes  the  PSF  based  upon  an  integration  of 
the  radial  spot  size  that  is  calculated  by  a  geometric  ray  trace.  If  the  aberrations  are 
twenty  times  the  diffraction  limit,  the  diffraction  method  will  independently  switch 
to  the  geometric  method. 
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IMAGE  SIMULATION:  DIFFRACTION  ABERRATIONS 

SIMULATING  THE  IMAGE  FORMED  BY  A  SINGLET  EYEPIECE 

OBJECT  HEIGHT  IS  13.8000  MILLIMETERS. 
ti.ee  nrt 
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CONFIGURATION  l  OF  1 
EXAMPLE  1.  H  SINGLE!  EYEPIECE, ZMX 

(a) 


(b) 


Figure  3.2.  (a)  A  PSF  grid  displaying  the  effects  of  distortion,  (b)  The  resulting  image 
simulation  by  convolving  an  image  source  file  with  the  PSF  grid. 


3.3  Zemax  Models 

Three  models  were  developed  in  this  thesis  for  all  the  scenarios  simulated.  A 
paraxial  model,  were  an  ideal  lens  refracts  all  rays  (for  all  wavelengths)  perfectly.  A 
chromatic  model,  were  multi-element  telephoto  lenses  are  modeled  to  demonstrate  the 
effects  of  chromatic  aberrations.  A  aberrated  model,  were  a  Zernike  phase  surface 
adds  controlled  amounts  of  aberrations.  These  models  incorporate  the  same  focal 
lengths  and  DVP  described  in  Section  2.8.1. 

3.3.1  The  Paraxial  Model. 

The  first  Zemax  model  created  uses  paraxial,  or  ideal,  lenses  for  the  development 
and  assessment  of  the  reconstruction  algorithm.  Paraxial  lenses  behave  ideally  and 
all  rays  converge  to  a  single  point  resulting  in  perfect  optical  system  response.  The 
model  is  shown  in  Figure  3.3. 

This  model  matches  the  theoretical  operation  as  described  in  literature,  thereby 
any  deficiencies  in  the  reconstruction  results  are  attributed  to  the  algorithm. 
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Figure  3.3.  Zemax  model  using  paraxial  lenses. 

3.3.2  Chromatic  Model. 

The  chromatic  model  uses  telephoto  lenses  modeled  after  prescriptions  from  patents [39, 
40] .  Due  to  the  difficulty  of  canvassing  through  thousands  of  patents,  these  prescrip¬ 
tions  are  not  an  exact  match  of  the  GCTEx  lenses.  The  number  of  lens  elements  differ 
and  the  F-numbers  are  larger;  however  the  focal  lengths  are  the  same.  The  goal  is 
to  have  something  representative  of  GCTEx  to  demonstrate  the  effects  of  chromatic 
aberrations.  Since  the  aperture  sizes  of  the  chromatic  model  lenses  are  smaller,  the 
focal  length  of  L2  was  increased  to  ensure  that  the  theoretical  Airy  disk  radius  was 
the  same  as  the  paraxial  model.  A  paraxial  lens  was  used  to  ensure  the  change  in 
longer  focal  lengths  did  not  negatively  impact  the  system  performance.  The  increased 
focal  length  of  Lo  lowers  the  angular  magnification  but  is  compensated  by  adjusting 
the  field  height  of  the  source  image  hie.  Figure  3.4  displays  the  chromatic  model. 

This  model  demonstrates  the  impact  of  chromatic  aberrations  on  the  reconstruc¬ 
tion  algorithm  and  the  wavelength  dependent  degradation  of  the  optical  system’s 
spatial  resolution. 
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Figure  3.4.  Zemax  model  using  telephoto  lenses. 

3.3.3  Zernike  Phase  Surface  Model. 

The  Zernike  Phase  Surface  function  in  Zemax  was  designed  to  model  optical  sys¬ 
tems  as  a  black  box  so  that  optical  designers  can  share  the  performance  of  the  system 
without  the  giving  the  exact  prescription.  The  surface  can  also  be  used  to  incorporate 
interferometric  measurements  of  optical  components  and  model  them  in  Zemax. 

The  model  uses  a  Zernike  Phase  Surface  to  control  the  amount  of  aberrations 
projected  onto  the  image  plane  from  the  focusing  lens.  The  Zernike  Phase  Surface 
modifies  the  incoming  wavefront  according  to  the  37  term  Zernike  FRINGE  polyno¬ 
mials.  However,  only  the  first  nine  Zernike  terms  are  used  to  represent  the  Siedel 
aberrations.  The  Zernike  Phase  Surface  is  placed  directly  after  a  paraxial  lens  which 
focuses  the  light  onto  the  detector.  The  telescope  portion  of  the  GCTEx  is  also 
made  of  paraxial  lenses  so  the  model  looks  exactly  like  the  paraxial  model  shown  in 
Figure  3.3. 
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3.4  Validation  of  Model  and  Simulation 


To  validate  the  paraxial  model  and  simulation,  the  radial  shift  given  by  Equa¬ 
tion  (2.11)  needs  to  match  the  deviation  angle  measured  by  Niederhauser.  A  monochro¬ 
matic  point  source,  modeled  as  a  perfect  circle,  centered  in  the  field-of-view  was  sim¬ 
ulated.  The  point  source’s  wavelengths  was  varied  from  400  nm  to  750  nm  in  25  nm 
increments.  The  deviation  angle  was  calculated  by  simulating  50  evenly-spaced  pro¬ 
jections  for  each  wavelength  increment.  Using  a  modified  script  from  Niederhauser, 
the  deviation  angle  is  computed  by  calculating  the  centroid  of  each  projection  and 
finding  a  circular  fit  for  all  of  the  centroids  as  shown  in  Figure  3.5. 


This  method  is  representative  of  how  the  measured  deviation  angles  were  obtained 
by  Niederhauser.  The  results  of  the  validation  are  shown  in  Figure  3.6.  As  expected, 
the  radial  shift  of  the  simulation  matches  the  theoretical  results. 
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Figure  3.6.  Measured  deviation  angle  from  a  simulated  point  source  at  multiple  wave¬ 
lengths. 


The  chromatic  model  is  validated  by  simulation  of  mercury  and  neon  pen  lamp 
point  sources  and  comparing  them  to  the  measured  lab  experiments.  The  chromatic 
simulation  uses  the  normalized  relative  intensities  from  Table  2.6  for  each  pen  lamp. 
The  model  was  focused  to  the  undeviated  wavelength  of  547.5nm,  near  the  center  of 
the  spectrum  of  interest.  The  size  of  the  simulated  spot  size  was  adjusted  to  match 
that  of  the  measured  data.  Figures  3.7(a)  and  3.7(b)  are  the  measured  and  simulated 
outputs  for  the  mercury  pen  lamp. 

Qualitatively,  the  model  matches  the  measured  data  well.  The  positions  of  each 
emission  line  match  with  each  other  and  the  blurring  on  the  435.8  nm  line  is  observed 
on  both.  The  404.6  nm  line  has  a  faint  signal  that  is  close  to  the  noise  floor,  making 
it  difficult  to  observe  on  the  measured  data  set.  The  lower  intensity  of  the  line,  the 
low  QE  response  at  that  wavelength,  detector  noise,  and  the  blurring  of  the  spot  are 
the  primary  causes  why  the  line  has  a  faint  signal. 
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(a) 


(b) 


Figure  3.7.  Measured  and  simulated,  respectively,  mercury  pen  lamp  data.  The  longer 
wavelengths  are  near  the  bottom  of  the  images. 

Using  the  same  configuration  as  the  mercury  pen  lamp,  the  neon  pen  lamp  was 
simulated  for  its  respected  wavelengths.  Figures  3.8(a)  and  3.8(b)  are  a  comparison 
of  the  measured  data  verse  the  simulated. 

The  simulation  also  models  the  GCTEx  chromatic  aberrations  well  for  the  longer 
wavelengths.  The  blurring  of  the  model  is  more  pronounced  at  the  longer  wavelengths; 
however,  considering  the  telephotos  are  not  an  exact  match  the  results  are  acceptable. 

The  Zernike  model  is  not  validated  experimentally  but  the  model  is  able  to  add 
controlled  amounts  of  aberrations  as  shown  in  Figures  3.9(a)-  3.9(d).  With  respect 
to  550nm,  10  waves  are  placed  on  each  of  the  Zernike  coefficients,  Z6  and  Z7,  to 
demonstrate  the  effects  of  coma. 

The  simulation  consists  of  four  projections  and  wavelengths  varying  from  400-750 
nm  in  50  nm  increments.  The  shorter  wavelengths  are  the  outermost  and  have  more 
pronounced  effects  from  the  aberrations  than  those  wavelengths  in  the  center. 
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(a) 
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Figure  3.8.  Measured  and  simulated,  respectively,  neon  pen  lamp  data.  The  longer 
wavelengths  are  near  the  bottom  of  the  images. 

3.5  Source  of  Transverse  Offset 

Determining  the  source  of  the  transverse  offset  is  important  so  that  it  can  be 
corrected  or  avoided  in  future  designs.  To  find  the  root  cause,  tilt  and  centration 
errors  for  each  surface  were  modeled.  Additionally,  the  Zernike  model  was  used  to 
determine  if  any  aberrations  may  be  the  cause.  The  main  criteria  for  finding  the  root 
was  looking  for  a  tilt  or  offset  that  caused  a  shift  in  the  undeviated  wavelength  when 
the  prism  is  rotated. 

No  tilts,  offsets,  and  aberrations  on  the  lenses  result  in  the  undeviated  wavelength 
to  shift  with  prism  rotation.  This  lead  to  investigating  the  prism  and  was  quickly 
identified  as  the  likely  candidate.  Hawks [21]  identified  the  transverse  offset  is  likely 
an  alignment  error  from  rotating,  or  “twisting”,  one  of  the  glasses  about  the  normal 
to  the  interface  surface  where  the  glasses  meet.  Figure  3.10  illustrates  an  exaggerated 
rotation  of  a  simple  top-down  or  bottom-up  perspective  of  the  prism,  separated  into 
the  two  glass  components  (red  and  blue). 

The  resulting  “twist”  causes  the  lateral  shift  and  phase  lag  observed  in  Figure  2.15. 
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Figure  3.9.  Four  orthogonal  projections  created  by  the  Zernike  model  displaying  the 
effects  of  coma.  A  uniform  spectral  intensity  point  source  sampled  at  25  nm  increments 
is  shown  for  each  projection.  The  “tail”  of  coma  is  more  pronounced  for  the  shorter 
wavelengths  (closer  to  the  edges  of  the  figures). 


52 


Figure  3.10.  A  top-down  or  bottom- up  perspective  of  an  exaggerated  rotation  about  the 
normal  to  interface  surface  where  the  glasses  (shown  by  the  red  and  blue  rectangles). 
This  rotation  is  believed  to  be  the  source  of  the  transverse  offset. 

Figure  3.11  is  a  Zemax  generated  footprint  diagram  of  orthogonal  projections  us¬ 
ing  the  mercury  pen  lamp  emission  lines.  The  546.8  nm  projections,  designated  by 
squares,  share  the  same  symmetric  offset  and  the  lag  in  phase  with  the  measured 
data. 

Figure  3.12  plots  the  transverse  offset  as  a  function  of  the  amount  of  prism  mis¬ 
alignment  using  the  Zemax  model.  A  85  mm  lens  and  a  detector  with  a  20  /im  pixel 
pitch  were  used  produce  the  amount  of  error  in  pixels.  The  transverse  error,  is  the 
lateral  shift  perpendicular  to  the  dispersion  path  while  the  wavelength  error  is  along 
with  the  path. 

The  results  of  this  simulation  are  used  to  provide  an  estimate  of  how  much  the 
prism  is  misaligned.  The  next  section  covers  how  the  transverse  offset  was  measured. 

3.5.1  Correcting  for  the  Transverse  Offset. 

The  first  step  towards  correcting  the  transverse  offset  was  to  determine  the  amount 
of  lateral  shift.  To  do  so,  a  elliptical-fit  MATLAB  script  developed  by  Niederhauser 
that  calculates  the  radius  and  eccentricity  of  elliptical  shapes  was  used.  From  mea- 
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Figure  3.11.  Footprint  diagram  showing  a  “twist”  in  the  prism  alignment.  Four  or¬ 
thogonal  projections  of  the  mercury  pen  lamp  wavelengths  are  shown.  The  546.8  nm 
projections,  designated  by  squares,  have  a  symmetric  offset  and  lag  in  phase  like  the 
measured  data  set  in  Figure  2.15.  [21] 
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Figure  3.12.  Transverse  and  wavelength  errors  for  varying  degrees  of  prism  alignment 
error.  Transverse  error  is  the  perpendicular  shift  from  the  dispersion  path  while  the 
wavelength  error  is  along  the  path.  [21] 
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sured  the  same  mercury  pen  lamp  data  used  to  observe  the  chromatic  aberrations, 
the  centroid  of  the  546.8  nm  line  from  each  projection  was  calculated  and  plotted  as 
shown  in  Figure  3.13.  These  centroids  are  then  input  into  the  elliptical-fit  script  for 
calculation. 
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Figure  3.13.  Centroids  from  each  projection  of  the  546.8  nm  line  of  a  mercury  pen  lamp 
were  used  to  determine  the  amount  of  lateral  shift.  The  elliptical-fit  script  calculates 
the  center,  radius,  and  its  eccentricity  of  the  shape.  The  radius  of  the  elliptical  shape 
must  be  accounted  for  in  the  reconstruction. 


The  results  from  the  elliptical-fit  script  computed  the  transverse  offset  as  16.9  ± 
0.4  pixels  and  was  nearly  a  perfect  circle.  From  Figure  3.12,  an  alignment  error  of 
20.5  arc  minutes  produces  a  transverse  offset  of  16  pixels.  A  20.5  arc  min  error  is 
reasonable  if  the  “twist”  tolerance  was  not  specified  during  fabrication  of  the  prism. 

To  correct  for  the  offset,  an  additional  shift  must  be  applied  to  account  for  the 
transverse  offset.  Figure  3.14  depicts  the  two  shifts  necessary  for  reconstruction,  the 
first  shift  is  the  uncorrected  shift  defined  in  Equation  (2.17)  which  places  the  object 
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onto  the  circumference  of  the  transverse  offset  and  the  second  shift  then  moves  the 
object  into  the  center  of  the  reconstructed  image. 


Figure  3.14.  The  first  shift,  described  by  the  STF,  places  the  object  on  the  circumfer¬ 
ence  of  the  transverse  offset.  A  second  shift  then  places  the  object  to  the  center  of  the 
reconstructed  image.  Each  shift  is  broken  into  its  corresponding  spatial  components. 

Applying  the  second  shift  to  each  element  of  the  STF,  Equation  (2.17)  now  be¬ 
comes 

Am)n(0  —  exp[  2?r£  •  Pshift\-  (3-1) 

The  projection  angle  shift,  p shift  is 

Pshift  =  {n  cos(2Trm/M)+Rsm(2'Km/M+(j)0ffset),n  sm(2Trm/M)—Rcos(2nm/M+(f)0ffset)), 

(3.2) 

where  R  is  the  radius  of  the  transverse  offset  and  <p0ffset  is  the  observed  phase  term 
for  each  wavelength  of  interest.  The  phase  term  is  found  empirically  using  calibrated 
light  sources.  Failure  to  account  for  the  transverse  offset  causes  the  shifted  projection 
not  to  align  properly  and  the  quality  of  the  reconstructed  image  degrades  significantly. 
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3.6  Summary 


In  this  chapter,  the  development  of  the  simulation  was  covered.  MATLAB  con¬ 
trols  functions  within  Zernax  to  produce  simulations  of  object’s  spatial  and  spectral 
content.  Three  Zemax  models  were  created:  the  paraxial,  chromatic,  and  zernike 
phase  surface.  Each  of  these  models  have  their  utility  in  assessing  different  aspects 
of  the  reconstruction  algorithm.  The  paraixal  model  was  validated  by  matching  the 
deviation  angles  measured  by  the  GCTEx  instrument.  The  chromatic  model  matches 
the  blurring  observed  in  the  GCTEx.  The  Zernike  phase  surface  demonstrated  that 
controlled  amounts  of  aberrations  can  be  injected  into  the  system.  This  chapter 
covered  the  source  of  the  transverse  offset  had  a  lateral  shift  of  16.9  pixels,  corre¬ 
sponding  to  a  20.5  arc  minute  alignment  error  within  the  prism.  Finally,  to  correct 
for  the  transverse  offset,  a  second  shift  in  the  reconstruction  algorithm  was  added. 
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4.  Reconstruction  Assessment 


With  the  models  and  simulations  set,  an  evaluation  of  the  reconstruction  algorithm 
is  required.  For  this  thesis,  only  the  fast  reconstruction  algorithm  by  Denting  [10] 
is  evaluated.  There  are  four  sets  of  test  cases  in  this  chapter,  a  set  for  each  of  the 
models  described  and  an  additional  set  taking  data  from  the  GCTEx  instrument. 

4.1  Metrics 

The  metrics  for  evaluating  reconstruction  performance  are  based  off  image  quality 
measures  used  to  gauge  performance  in  compressing  digital  images.  Since  compression 
and  the  reconstruction  process  lose  information  original  content,  the  image  quality 
metrics  are  suited  for  measuring  the  reconstruction’s  performance.  There  are  two  cat¬ 
egories  of  image  quality  metrics  used:  pixel-difference-based  and  edge  based  measures 
[1,  8], 

A  common  pixel-difference-based  metric  is  the  Mean  Square  Error  (MSE).  For  an 
image  where  there  are  M  and  N  number  of  pixels  in  the  row  and  column  directions, 
respectively,  the  MSE  is  given  by 

M  N 

MSE  ~  MN  ££(^-y,)2,  (4.D 

j= 1  k= 1 

where  Xj^  denotes  pixel  (j,  k)  of  the  original  image  and  x']  k  pixel  (j,  k )  of  the  recon¬ 
structed  image.  Higher  values  of  MSE  corresponds  to  poor  image  quality. 

The  Laplacian  Mean  Square  Error  (LMSE)  has  emphasis  on  the  edges  of  the  scene. 
This  metric  will  gage  the  sharpness  of  the  edges,  offset  in  edge  position,  and  missing 
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edges.  The  LMSE  is  defined  as 


LSME 


Ef.iEhl  [0(Hk)-0(x'j:t) ]■ 


EiELOfe 


N 


(4.2) 


where  the  Laplacian  operator,  0(xjtk)  is 


O (•£ j,k)  *£(j+i),fc  T  ^(j— i),fc  T  1)  T  •^j,(k— i)  dx'j  /,, .  (4.3) 

A  larger  value  corresponds  to  a  poor  image  quality.  Together,  these  metrics  will 
provide  an  objective  method  of  measuring  blurring  or  signal  loss  in  the  reconstructed 
images. 

4.2  Paraxial  Model  Test  Cases 

The  ideal  optical  system  test  cases  use  the  paraxial  model  to  isolate  the  recon¬ 
struction  algorithm  as  the  source  for  any  degradation  in  image  quality.  Simulations 
evolve  from  a  simple  scene  illuminated  by  a  mercury  pen  lamp,  to  a  static  scene,  and 
then  a  polychromatic  scene. 

4.2.1  Parameters  for  Reconstruction. 

A  simple  scene  with  distinct  spectral  features  is  the  first  test  case  scenario.  A  5x5 
grid  of  the  alphabet  is  illuminated  by  a  mercury  pen  lamp.  The  mercury  pen  lamp 
produces  spectrally  narrow  emission  lines,  therefore  it  is  assumed  the  width  of  the 
emission  lines  are  negligible.  Four  of  the  five  strong  mercury  emission  lines  in  the 
visible  region  are  used,  the  579.0  nm  line  is  ignored  due  to  the  close  proximity  of  the 
576.9  nm  line.  Figure  4.1  is  the  result  of  one  projection  simulated.  Note  that  each 
wavelength  has  been  scaled  according  to  Table  2.6. 

The  fast  reconstruction  algorithm  has  a  filtering  parameter,  /i,  which  provides 
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Figure  4.1.  A  5x5  grid  of  letters  illuminated  by  a  mercury  pen  lamp.  The  wavelengths 
are  scaled  according  to  their  relative  intensities. 


stability.  Different  values  of  the  filtering  parameter  were  tested,  ultimately  leading 
to  the  conclusion  that  //  must  be  greater  than  or  equal  to  one  for  stability.  Setting 
/i  much  larger  than  one  essentially  applies  no  filter  and  the  results  are  similar  to  the 
unfiltered  response.  The  filtering  also  forces  the  reconstructed  images  to  have  positive 
and  negative  values,  centered  about  zero.  The  user  must  select  which  range  provides 
the  best  results,  and  then  set  the  other  half  equal  to  zero.  If  the  negative  range  was 
selected,  the  absolute  value  was  taken  to  provide  more  physical  results. 

The  results  of  reconstructing  the  mercury  pen  lamp  are  shown  in  Figures  4.2(a)- 
4.2(d).  For  each  wavelength  of  the  mercury  pen  lamp,  the  unfiltered  (left)  and  the 
filtered  (right)  reconstructed  image  is  shown.  This  first  set  of  results  are  from  100 
projections  simulated.  The  404  nm  and  576  nm,  Figures  4.2(a)  and  4.2(d)  respectively, 
the  negative  range  of  the  filtered  reconstructed  image  were  chosen. 

As  shown  in  Figures  4.2(a)  and  4.2(d),  the  reconstruction  algorithm  is  able  to 
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Figure  4.2.  (a-d)  Displays  the  reconstructed  image  for  the  unfiltered  (left)  and  filtered 
(right)  of  each  wavelength  of  the  mercury  pen  lamp. 
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piece  together  the  spectral  contents  well.  The  filtered  images  clean  the  surrounding 
noise  but  at  the  cost  of  signal  drop.  To  calculate  the  signal,  the  reconstructed  images 
must  be  scaled  by  the  number  of  projections  used  in  reconstruction.  Since  half  of  the 
filtered  range  is  chopped,  it  must  be  multiplied  by  two  to  achieve  signal  strengths  on 
par  with  the  unfiltered  images.  A  different  multiplication  factor  is  required  if  using 
values  outside  of  [i  —  1.  Figure  4.3  plots  the  mean  signal  of  the  grid  (not  the  entire 
image)  for  each  wavelength  and  compares  it  with  input  signal.  The  mean  of  the  grid 
is  computed  because  the  reconstructed  images  do  not  have  a  uniform  response  across 
the  grid. 


Figure  4.3.  The  mean  signal  of  the  grid  for  the  filtered  and  unfiltered  reconstructed 
images.  The  filtered  images  are  have  a  large  drop  where  the  input  signal  is  strong. 


Where  the  input  signal  is  at  its  peaks,  the  filtered  signals  are  significantly  lower, 
up  to  60%  lower  than  the  input  signal.  The  MSE  is  used  to  quantify  the  quality  of 
the  images  are  shown  in  Figure  4.4. 
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Figure  4.4.  MSE  for  the  unfiltered  and  filtered  reconstructed  images.  With  the  excep¬ 
tion  of  the  435  nm  line,  the  filtered  response  has  a  better  image  quality. 

Figure  4.4  shows  the  filtered  images  have  higher  image  quality  with  the  exception 
of  the  435  nm  line.  As  mentioned  earlier,  the  filtered  435  nm  reconstructed  mean 
signal  is  substantially  lower  than  the  input  signal  and  as  a  result  the  image  quality 
is  degraded  as  well. 

The  image  quality  is  dependent  on  the  number  of  projections  used  in  the  recon¬ 
struction  algorithm.  Figure  4.5  plots  the  MSE  for  un filtered  images  with  the  number 
of  projections,  evenly  spaced  from  0  to  27 r,  are  increased. 

This  shows  that  a  minimum  of  64  projections  are  required  to  obtain  the  best 
reconstructed  images.  After  64  projections,  the  marginal  improvement  is  effectively 
zero  for  all  the  wavelengths. 
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Figure  4.5.  The  performance  of  the  reconstruction  varies  with  the  number  of  projec¬ 
tions  used.  The  reconstruction  performance  will  not  increase  after  64  projections. 


4.2.2  Static  Scene. 

A  more  complex  scene  is  used  to  determine  if  the  reconstruction  algorithm  can 
extract  spectral  features  of  a  scene.  A  hypercube  taken  from  the  AVIRIS  platform 
database  [26]  was  simulated.  This  scenario  demonstrates  the  effects  of  the  chromo- 
tomographic  process  and  reconstruction  algorithm  data  taken  from  a  calibrated  HSI 
instrument. 

The  AVIRIS  spectral  range  covers  224  bins  with  a  resolution  of  approximatelylO 
nm  for  each  bin.  AVIRIS  frames  from  the  400-740  nm  range  were  selected.  For  each 
bin,  a  256x256  pixel  clip  of  the  scene  was  selected  and  zero  padded  to  a  768x768  hie 
size.  This  zero-padding  is  necessary  to  ensure  the  256x256  scene  was  not  truncated 
from  the  dispersion  of  the  prism.  The  image  held  height  was  set  such  the  output 
simulations  matched  the  same  dimensions  as  the  inputted  image  so  it  is  easily  com¬ 
parable.  Each  frame  was  normalized  by  the  maximum  bit-size  of  the  AVIRIS  format. 
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A  projection  simulated  at  0  degrees  is  shown  in  Figure  4.6. 


Figure  4.6.  Simulation  of  the  AVIRIS  scene  at  a  single  projection.  A  zero-padding  was 
adding  to  ensure  all  information  was  collected. 


A  total  of  36  projections  were  simulated  and  placed  through  the  reconstruction 
algorithm.  A  high-pass  filter  was  added  to  both  the  unfiltered  and  filtered  reconstruc¬ 
tion  to  remove  the  lower  spatial  frequency  content  that  was  degrading  the  contrast 
of  the  scene.  Figures  4.7(a)-4.7(c)  shows  the  results  of  the  reconstruction  at  the 
beginning,  middle,  and  end  of  the  band. 

The  unfiltered  reconstruction  clearly  outperforms  the  filtered  method  at  the  mid¬ 
dle  wavelengths.  However,  the  filtered  reconstruction  performs  better  than  the  longer 
wavelengths.  Considering  the  entire  data  set,  the  longer  wavelengths  produced  better 
quality  reconstructed  images.  Shorter  wavelengths  lack  of  performance  are  attributed 
to  the  cone  of  missing  information,  i.e.  the  volume  of  the  cone  is  larger  for  the  shorter 
wavelengths.  The  middle  wavelengths  are  inundated  with  the  all  of  the  other  spectral 
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Figure  4.7.  (a-c)  Displays  the  reconstructed  image  for  the  unfiltered  (left)  and  filtered 
(right).  The  center  image  is  from  the  AVIRIS  hypercube  at  the  given  wavelength. 
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Figure  4.8.  Demonstrates  that  the  reconstruction  process  can  identity  spectral  features 
such  as  when  water  begins  to  absorb  light  at  700  nm.  Three  successive  bins  at  (a)  695 
nm,  (b)  705  nm,  and  (c)  715  nm  are  shown.  The  image  on  the  left  is  from  the  unfiltered 
method,  the  center  image  is  from  the  AVIRIS  hypercube,  and  the  filtered  method  is 
on  the  right. 


content,  its  reconstructed  image  is  thereby  degraded. 

Both  the  filtered  and  unfiltered  reconstruction  algorithms  are  able  to  pick  out 
the  distinct  spectral  features  such  as  the  water  absorption  line.  Figures  4.8(a)-4.8(c) 
show  three  bands  that  show  when  water  starts  to  absorb  at  700  nm. 

The  dark  patches  in  the  scene  are  vegetation  and  as  water  starts  to  absorb  the  veg¬ 
etation’s  reflectance  is  increased  thereby  the  reflected  signal  from  vegetation  increases. 
These  test  demonstrated  the  fast  reconstruction  algorithm  is  able  to  reconstruct  the 
longer  wavelengths  adequately;  however,  more  advanced  reconstruction  algorithms 
are  necessary  to  fill  in  the  cone  of  missing  information. 

4.2.3  Polychromatic  Scene. 

A  polychromatic  scene  is  modeled  by  four  Gaussian  sources,  representative  of  a 
light-emitting-diode,  at  varying  wavelengths.  Each  of  these  sources  are  then  spatially 
tied  to  a  letter  in  the  acronym,  AFIT,  as  shown  in  Figure  4.9.  The  letter  A  is 
associated  with  the  red  source,  letter  F  with  the  blue,  letter  1  with  green,  and  the 
letter  T  with  yellow. 

The  sources  are  arranged  such  that  they  overlap  spectrally  but  not  spatially.  Each 
source  has  a  maximum  intensity  of  one  and  a  FWHM  of  50  nm.  The  blue,  green, 
and  yellow  sources  intersect  at  their  corresponding  FWHM,  while  the  red  source  is 
separated  further  to  cover  more  of  the  visible  spectrum.  Figure  4.10  plots  the  spectral 
distribution  of  for  all  the  sources. 

Like  the  static  scene,  the  size  of  the  image  hie  is  zero-padded  to  ensure  that  all  of 
the  dispersed  scene  is  captured  and  rescaled  for  comparison.  The  granularity  of  the 
spectral  profile  is  in  10  nm  increments  and  a  total  of  36  projections  were  simulated. 
Figure  4.11  shows  the  simulation  output  from  a  projection  at  0  degrees.  Note  that 
the  letters  are  dispersed  at  different  distances  onto  the  FPA. 


Figure  4.9.  Illustration  of  the  scene  simulated,  with  no  spatial  overlap  but  spectral 
overlap  from  each  source. 


Figure  4.10.  Spectral  profile  of  the  simulation.  The  FWHM  of  each  Gaussian  source  is 
50  nm  and  spaced  to  overlap  one  another. 
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Figure  4.11.  Each  letter  is  shifted  at  different  distances  due  to  their  varying  spectral 
content.  The  latter  three  letters  have  approximately  50  nm  of  overlap. 

The  reconstructed  images  at  their  respected  peaks  are  shown  in  Figures  4.12(a)- 
4.12(d).  For  each  figure,  the  unfiltered  and  filtered  reconstructed  images  are  shown 
left  and  right,  respectively.  As  depicted,  the  filtering  reduces  the  noise  from  the 
surrounding  spectral  content. 

The  spectrum  of  the  reconstructed  images  is  computed  as  the  mean  intensity  for 
each  spatial  letter  as  a  function  of  wavelength.  Figures  4.13(a)  and  4.13(b)  show 
the  peaks  at  the  shorter  wavelengths  are  shifted  and  that  the  intensity  increases 
considerably  at  longer  wavelengths. 

The  unfiltered  spectrum  has  significant  bleeding  into  the  blue  source  from  the 
other  spectral  channels.  The  filtered  spectrum  considerably  eliminates  most  of  the 
cross-source  bleeding.  The  unfiltered  spectrum  results  have  much  larger  FWHM  for 
each  source,  as  summarized  in  Table  4.1.  The  filtered  spectrum’s  FWHM  are  close 
to  the  actual  input  FWHMs  of  50  nm.  The  blue  and  green  source  peaks  are  shifted 
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(a)  450  nm 


(b)  500  nm 


(c)  550  nm 


(d)  650  nm 


Figure  4.12.  (a-d)  Displays  the  reconstructed  image  for  the  unfiltered  (left)  and  filtered 
(right)  of  each  spectral  peak  wavelength. 
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Figure  4.13.  (a)  Unfiltered  spectrum  showing  bleeding  of  the  sources,  (b)  Filtered 

spectrum  significantly  reduces  the  spectral  widening.  The  shorter  wavelength  peaks 
shift  as  the  intensity  of  the  longer  wavelengths  increases. 


by  10  nm  for  both  sets  of  reconstructed  images. 

All  of  the  sources  have  lower  reconstructed  intensities  than  the  actual  source 
amount.  The  drop  in  intensity  is  more  pronounced  at  the  shorter  wavelengths  where 
the  dispersion  from  the  prism  is  at  its  greatest.  The  fast  reconstruction  algorithm  is 
unable  to  maintain  the  absolute  radiometric  intensity  but  the  general  profile  or  shape 
of  the  spectrum  is  maintained. 


Table  4.1.  FWHM  of  the  unfiltered  and  filtered  reconstructed  spectrum.  Ideally,  the 
FWHM  should  be  50  nm  for  each  source. 


Letter  (color) 

Unfiltered  FWHM  (nm) 

Filtered  FWHM  (nm) 

A  (red) 

97.0 

65.3 

F(blue) 

66.4 

46.9 

I  (green) 

65.8 

52.0 

T  (yellow) 

73.7 

52.0 

The  MSE  of  the  filtered  images  are  also  better  than  the  unfiltered  as  indicated 
in  Figure  4.14.  Interestingly,  MSE  degrades  at  the  spectral  peaks,  more  so  for  the 
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shorter  wavelengths  than  the  longer  wavelengths. 


Figure  4.14.  MSE  spikes  at  the  shorter  wavelength  peaks  due  to  the  difference  of  the 
reconstructed  signal  strength  verse  the  actual  signal.  The  error  in  the  filtered  recon¬ 
structed  images  is  lower  because  the  residual  noise  artifacts  from  the  other  spectral 
channels  are  reduced. 

The  elevated  MSE  at  the  peaks  are  due  to  the  difference  from  the  reconstructed 
signal  strength  verse  the  actual  signal.  As  the  actual  signal  peak  occurs,  the  large 
signal  difference  drives  the  MSE  upwards.  Conversely,  as  the  wavelengths  move  from 
the  spectral  peaks,  the  difference  between  signals  is  not  as  great  resulting  in  a  lower 
MSE.  This  explains  why  the  MSE  at  the  650  nm,  where  the  reconstructed  intensity 
is  greater,  has  a  lower  MSE  value. 

In  short,  the  fast  reconstruction  algorithm  is  able  to  resolve  the  spectral  features; 
however,  it  is  unable  to  maintain  its  absolute  radiometric  quantity.  Filtering  is  re¬ 
quired  to  prevent  noise  from  the  other  spectral  channels  from  seeping  into  the  channel 
of  interest. 
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4.3  Chromatic  Model  Test  Case 


This  test  case  uses  the  chromatic  model  to  investigate  how  chromatic  blurring 
impacts  reconstruction.  A  simple  scene  with  an  object  with  high  spatial  frequency, 
shown  in  Figure  4.15,  is  simulated. 


Figure  4.15.  This  object  has  high  spatial  frequency  which  will  be  blurred  by  chromatic 
aberrations. 

The  object’s  spectral  intensity  is  uniform  over  the  spectral  range  from  400-750 
nm  and  the  scene  is  simulated  in  25  nm  increments.  The  chromatic  model  was 
focused  at  the  550  nm  wavelength.  A  uniform  spectrum  was  chosen  to  remove  any 
negative  impacts  the  spectral  features  may  have  in  the  reconstruction  algorithm. 
Since  high  spatial  frequency  content  is  lost  due  the  cone  of  missing  information, 
the  same  test  case  was  also  simulated  in  the  paraxial  model  for  comparison.  The 
filtered  reconstruction  does  not  perform  well  reconstructing  the  high  spatial  frequency 
content,  so  only  the  unfiltered  results  are  presented.  Figures  4.16(a)-4.16(d)  show  the 
reconstructed  images  for  paraxial  (left)  and  chromatic  (right)  simulations  for  several 
wavelengths  over  the  spectral  range. 
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(c)  650nm 


(d)  750nm 


Figure  4.16.  (a-d)  Reconstruction  of  the  high  spatial  frequency  object  using  the  parax¬ 
ial,  left,  and  chromatic,  right,  models.  Due  to  the  cone  of  missing  information,  both 
sets  of  reconstructed  images  have  a  significant  degradation  in  image  quality.  Loss  in 
image  quality  is  more  pronounced  at  the  ends  of  the  spectral  range. 
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The  400  nm  reconstructed  images  do  not  constructively  combine  each  shifted 
projection  as  effectively  as  the  other  wavelengths,  resulting  in  the  residual  artifacts 
to  be  more  prevalent.  This  is  consistent  with  the  static  scene  and  AFIT  test  cases 
where  the  short  wavelength  reconstructed  images  have  poor  image  quality  and  the 
longer  wavelengths  performed  better.  This  again  is  attributed  to  the  cone  of  missing 
information,  where  the  longer  dispersion  at  the  short  wavelengths  cause  a  larger  cone 
to  be  formed  therefore  degraded  reconstructed  images. 

The  chromatic  aberrations  shift  the  spectral  location  of  the  highest  quality  recon¬ 
structed  image  as  shown  in  Figures  4.16(b)  and  4.16(c).  Since  the  chromatic  model 
was  focused  at  550  nm,  the  highest  quality  reconstructed  image  occurs  near  the  550 
nm  while  the  paraxial  model  forms  the  best  reconstructed  image  at  650  nm.  To  qual¬ 
ify  the  image  quality,  the  LMSE  metric  is  used  to  put  emphasis  edges  of  the  high 
spatial  frequency  object.  Figure  4.17  plots  the  LMSE  over  the  spectral  range. 

The  chromatic  model  performs  better  than  the  paraxial  model  at  the  shorter 
wavelengths  likely  from  the  adjustment  in  focus  in  the  chromatic  model.  However, 
the  chromatic  model  performs  poorly  at  the  longer  wavelengths  and  a  the  chromatic 
model  substantial  difference  in  the  troughs  of  the  curves.  The  chromatic  model  has  a 
narrow  operating  range  from  roughly  525  nm  to  575  nm  where  high  spatial  frequency 
content  is  best  maintained. 

4.4  Zernike  Phase  Surface  Model  Test  Case 

This  test  case  is  similar  to  the  chromatic  model  with  the  exception  that  different 
combinations  of  Zernike  coefficients  are  used  to  add  controlled  amounts  of  aberra¬ 
tions.  Like  the  chromatic  model  test,  each  aberration  term  is  simulated  by  a  uniform 
spectrum  to  observe  how  the  aberrations  degrade  image  quality  of  the  high  spatial 
frequency  object,  shown  in  Figure  4.15.  This  provides  a  characterization  of  which 
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Figure  4.17.  The  high  spatial  frequency  content  is  best  preserved  at  550  nm  where  the 
focus  of  the  chromatic  model  was  at.  A  small  operating  range  exists  for  maintaining 
the  spatial  content. 


aberrations  are  the  most  detrimental  for  the  future  spaced-based  system  to  avoid. 

This  test  case  only  considers  the  effects  of  astigmatism,  coma,  and  spherical  aber¬ 
rations.  From  Table  2.3,  the  astigmatism  and  coma  aberration  terms  consist  of  two 
Zernike  coefficients,  each  of  the  coefficients  are  equally  increased.  Furthermore,  the 
coefficients  where  balanced  such  that  only  the  aberration  under  test  existed  in  the 
Zernike  phase  surface.  For  example,  coma  consists  of  the  Z§  and  Z-j  coefficients  which 
are  also  present  in  the  tilt  combination.  To  zero  out  the  tilt,  an  offsetting  amount  of 
the  Z\  and  Z2  were  applied. 

Each  aberration  was  incrementally  increased  from  1  wave  to  5  waves  to  10  waves. 
All  of  the  aberration  terms  are  referenced  to  550  nm.  Zemax’s  geometric  image 
simulation,  vice  the  image  simulation,  was  used  for  these  simulations  since  geometric 
ray  tracing  is  required  in  the  presence  of  strong  aberrations.  Once  again,  the  LMSE 
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metric  is  used  to  qualify  the  high  spatial  frequency  object.  Figures  4.18(a)-  4.18(c) 
plot  the  LMSE  for  each  of  the  aberrations  as  the  number  of  waves  increases.  Each 
plot  contains  a  reference  LMSE  where  no  aberrations  were  added. 


(a)  1  wave  (b)  5  waves 


(c)  10  waves 


Figure  4.18.  (a-c)  LMSE  for  each  aberration  as  the  number  of  waves  are  increased. 
One  wave  of  aberrations  have  little  impact  on  the  image  quality.  However,  steadily 
increasing  the  amount  of  waves  shows  that  coma  degrades  the  reconstruction  perfor¬ 
mance  the  most.  The  longer  wavelengths  reconstruction  performance  degrades  more 
rapidly  than  the  shorter  wavelengths. 


A  subjective  LMSE  threshold  of  0.990  is  used  determining  if  the  reconstruction 
images  are  acceptable  for  identifying  the  actual  object.  The  threshold  was  based  upon 
looking  at  reconstructed  images  from  the  simulations  produced  without  any  aberra¬ 
tions.  Spatial  features,  such  as  the  spokes  of  the  object,  are  clearly  discernible  with 
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a  LMSE  of  0.990  or  lower.  Figures  4.19(a)  and  4.19(b)  illustrate  two  reconstructed 
images  across  the  LMSE  threshold. 


Ill  !  II  III  III  ill  III  HI  III  III  1 1 1 1  III  III  III  III  ill  ill  HI  III  III  1 1 1 1 


(a)  LMSE=0.9916  (b)  LMSE=0.9873 

Figure  4.19.  Two  reconstructed  images  above  (a)  and  below  (b)  the  subjective  LMSE 
threshold  value  of  0.990. 

Once  again,  the  shorter  wavelengths  perform  worse  than  the  longer  wavelengths. 
The  550  nm  reconstructed  image  produces  the  best  reconstructed  image  for  all  sets 
of  reconstructed  images.  This  differs  from  the  chromatic  test  case  where  the  paraxial 
model  produced  the  best  reconstructed  at  650  nm.  The  results  vary  from  differences 
in  the  image  simulations  as  the  geometric  image  simulation  does  not  maintain  the  ob¬ 
ject’s  high  spatial  frequency  content  as  well  as  the  convolution  image  simulation.  The 
differences  in  the  simulations  are  trivial  provided  trends  of  how  aberrations  impact 
reconstruction  performance  are  identified. 

One  wave  from  each  of  the  three  aberrations  have  little  effect  on  the  reconstruc¬ 
tion  performance.  However  increasing  the  amount  to  5  waves,  the  reconstruction 
performance  is  degraded  predominantly  at  the  longer  wavelengths.  Coma  clearly  neg¬ 
atively  effects  the  reconstruction  performance  the  most,  while  spherical  aberrations 
degrade  image  quality  slightly  more  than  astigmatism.  At  10  waves  of  aberrations, 
the  negative  effects  on  reconstruction  are  compounded.  Coma  effectively  blurs  the 
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reconstructed  images  over  entire  spectral  range.  Astigmatism  and  spherical  aberra¬ 
tions  have  a  narrow  range  of  40  nm  centered  at  550  nm  where  the  image  quality  is 
below  the  threshold  value. 

4.5  Lab  Experiment  Test  Cases 

Simple  laboratory  experiments  of  a  monochromatic  scene  of  the  Air  Force  bar  chart 
was  collected.  Next,  an  iris  directly  in  front  of  a  mercury  pen  lamp  was  placed  to 
simulate  a  point  source.  For  each  of  these  experiments,  the  transverse  offset  correction 
was  applied  to  the  reconstruction  algorithm. 

4.5.1  Air  Force  Bar  Chart. 

An  Air  Force  bar  chart  was  measured  using  a  mercury  light  source  with  a  narrow- 
band  filter  at  546.8  nm.  The  prism  was  spinning  at  a  rate  less  than  1  rev/sec.  Data 
recorded  at  a  single  projection  is  shown  in  Figure  4.20.  Notice  that  two  bar  charts  are 
observed  in  the  measured  projection,  with  one  in  focus  and  the  other  blurred.  The 
blurred  bar  chart  is  believed  to  be  a  ghost  reflection  caused  by  internal  reflections 
within  the  instrument. 

It  is  possible  that  secondary  blurred  bar  chart  is  from  another  emission  line  from 
the  source;  however,  the  shift  between  each  bar  chart  is  more  than  90  pixels  which 
corresponds  to  1.2  degrees  of  deviation.  Since  the  bar  chart  in  focus  is  at  546.8  nm, 
the  1.2  degrees  of  deviation  corresponds  to  emission  lines  roughly  at  460  nm  or  800 
nm,  neither  of  which  are  strong  emission  lines  for  a  mercury  source. 

A  total  of  50  projections  of  the  bar  chart  were  collected  and  applied  the  recon¬ 
struction  algorithm  without  any  correction  for  transverse  offset.  The  resulting  recon¬ 
structed  image  is  shown  in  Figure  4.21. 

The  reconstructed  algorithm  severely  degrades  the  image.  The  transverse  offset 
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Figure  4.20.  Air  Force  bar  chart  illuminated  by  a  monochromatic  source  at  546. 8nm. 


causes  the  object  to  destructively  interfere  with  itself  resulting  in  a  blurred  image. 
This  is  quickly  resolved  by  applying  change  in  the  reconstruction  algorithm,  described 
in  Section  3.5.1,  to  account  for  the  transverse  offset.  The  reconstructed  image  of  the 
bar  chart  is  shown  in  Figure  4.22. 

A  significant  improvement  from  the  two  reconstructed  images.  The  spatial  content 
of  the  image  is  now  maintained  with  the  correction.  The  lower  spatial  frequencies 
are  saturating  the  reconstructed  image  thereby  lowering  the  contrast.  A  Gaussian 
high-pass  filter  is  also  applied  to  the  reconstructed  image  to  enhance  the  contrast  as 
shown  in  Figure  4.23. 

The  MTF  of  the  corrected  reconstructed  image  with  and  without  the  high-pass 
filter  was  calculated  using  Equation  (2.5)  by  Hawks  [20].  The  results  are  plotted  in 


81 


Figure  4.21.  Air  Force  bar  chart  reconstructed  without  applying  the  transverse  offset 
correction. 

Figure  4.24. 

The  high  pass  filter  improves  the  MTF  by  4  times  the  lower  frequencies.  The 
significance  of  these  results  validates  that  transverse  offset  correction  and  upgrade  in 
electronics  deliver  the  correct  prism  angles  necessary  to  perform  accurate  reconstruc¬ 
tion.  No  additional  corrections  were  applied  or  manual  manipulation  of  the  data  was 
used  in  the  data  collection  and/or  reconstruction  algorithm. 

4.5.2  Mercury  Pen  Lamp. 

The  objective  of  this  test  case  is  to  determine  if  the  fast  reconstruction  algorithm 
produce  a  data  cube  of  a  spectrally  simplistic  scene.  A  point  source  using  a  mercury 
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Figure  4.22.  Transverse  offset  applied  and  the  bar  chart’s  lines  are  now  observable  and 
distinct. 


pen  lamp  was  collected  using  the  GCTEx  instrument.  The  reconstructed  images 
of  the  mercury  pen  lamp  were  constructed  from  50  projections  and  are  shown  in 
Figures  4.25(a)-  4.25(c).  The  instrument  was  adjusted  such  that  the  546.8  nm  line 
was  in  focus. 

The  transverse  offset  correction  was  applied  and  the  phase  lag  terms  were  adjusted 
so  each  spectral  bin  would  converge  to  a  point.  Table  4.2  lists  amount  of  phase  lag  for 
required  for  each  emission  line.  All  of  the  phase  lags  terms  are  calculated  empirically. 

The  435.8  nm  spot  size  is  considerably  larger  than  the  other  emission  lines  due  to 
chromatic  aberrations.  The  phase  lag  correction  is  not  perfect  as  the  reconstruction 
doubled  the  spot  size  of  the  435.8  nm  line.  However,  it  does  perform  well  where  the 
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Figure  4.23.  High  pass  filter  improves  the  contrast  of  the  reconstructed  image.  [20] 

emission  lines  had  a  tighter  focus.  The  results  show  the  corrected  algorithm  is  capable 
of  adjusting  for  scenes  with  varying  spectral  content. 

To  provide  an  estimate  of  the  fast  the  reconstruction  algorithm  is,  this  test  case 
was  clocked  using  MATLAB’s  code  timing  mechanisms.  To  load  the  512x512  pixels 
camera  hies,  read  in  the  prism  angle  data,  backproject  50  projections  for  3  spec¬ 
tral  bins,  and  apply  a  filter  takes  an  average  of  3.50  seconds.  The  previous  AFIT 
backprojection  algorithm  required  several  minutes  for  an  equivalent  reconstruction. 
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Figure  4.24.  MTF  calculated  from  the  US  Air  Force  bar  chart  with  and  without  the 
high-pass  filter.  The  high-pass  filter  substantially  improves  the  resolving  power  of  the 
reconstruction  algorithm.  [20] 


Table  4.2.  Phase  lag  adjustment  required  for  each  emission  line  in  reconstruction. 


Emission  Line  (nm) 

Phase  Lag  (Degrees) 

435.8 

-35.21 

546.8 

-11.92 

576.9 

0 

4.6  Summary  of  Results 

This  chapter  demonstrated  that  the  fast  reconstruction  algorithm  is  capable  of 
reconstructing  simple  and  complex  scenes.  For  each  reconstructed  image,  64  projec¬ 
tions  that  span  from  0  to  27r  produce  the  best  results.  More  than  64  projections  do 
not  provide  any  significant  improvement  in  the  reconstruction  results. 

The  fast  reconstruction  algorithm  is  able  to  reconstruct  static  scenes  and  distin¬ 
guish  spectral  features  such  as  the  water  absorption  line.  However,  the  algorithm’s 
performance  over  the  spectral  range  is  poor  due  to  the  cone  of  missing  informa- 
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(c)  576.9  nm 


Figure  4.25.  (a-c)  Reconstructed  images  of  the  mercury  pen  lamp  point  source 


tion.  The  reconstruction  algorithm  is  able  to  identify  the  shape  of  spectral  sources 
as  demonstrated  with  the  polychromatic  test  case.  However,  shifts  in  the  peak  along 
with  an  lowered  intensities,  particularly  at  the  shorter  wavelengths,  prevent  it  from 
providing  absolute  radiometric  profile.  The  fast  reconstruction  algorithm  has  utility 
as  a  “quick  look”  algorithm  and  can  an  input  for  more  complex  algorithms  such  as 
the  algorithm  developed  by  Gould  [16]. 

From  the  chromatic  test  case,  high  spatial  frequency  objects  do  not  maintain 
integrity  attributed  to  the  cone  of  missing  information.  The  reconstruction  perfor¬ 
mance  is  degraded  by  the  effects  of  chromatic  and  monochromatic  aberrations.  The 
chromatic  aberrations  force  a  small  operating  range  where  the  blurring  does  not  sig¬ 
nificantly  degrade  reconstruction  performance.  The  aberration  test  case  concluded 
that  coma  hindered  reconstruction  performance  the  most. 

Finally,  lab  experiments  showed  accounting  for  the  transverse  offset  is  critical 
for  the  system  to  obtain  valid  reconstructed  images.  The  bar  chart  test  validated 
that  the  upgrade  in  electronics  now  sufficiently  sync  the  camera  and  prism  angles. 
The  instrument  and  reconstruction  algorithm  are  ready  for  more  complex  scenes  and 
transients  to  be  observed. 
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5.  Conclusions  and  Recommendations 


This  chapter  provides  a  summary  of  the  research  conducted  in  this  thesis.  Conclu¬ 
sions  are  made  on  the  model  and  simulation,  the  fast  reconstruction  algorithm,  and 
the  transverse  offset.  Future  work  recommendations  are  then  given.  The  conclusions 
made  do  not  fully  summarize  the  instrument’s  performance.  Since  the  instrument’s 
performance  is  dependent  on  the  scene,  the  conclusions  drawn  are  specifically  for  the 
scenes  simulated. 

5.1  Model  and  Simulation 

Zemax  and  MATLAB  proved  to  be  a  powerful  combination.  The  simulated  scenar¬ 
ios  are  a  small  portion  of  the  simulation’s  capability.  The  simulated  scenarios  in  this 
thesis  demonstrated  the  orthogonal  basis  set  of  possible  scenarios:  spatially  complex 
but  spectrally  simple  and  spectrally  complex  but  spatially  simple.  Furthermore,  the 
simulation  allows  flexibility  in  the  design  of  the  optical  system.  Controlled  amounts  of 
aberrations  or  optical  alignment  errors  can  be  analyzed  for  future  trade-space  choices 
in  the  design  of  the  space-based  system. 

Evaluation  of  the  Zemax  models  demonstrated  that  the  realistic  models  of  the 
GCTEx  can  be  made.  The  paraxial  model  provides  a  tool  for  simulating  the  theoret¬ 
ical  limits  of  the  chromtomographic  instrument.  Additionally,  it  was  an  invaluable 
tool  for  developing  the  reconstruction  algorithm  and  identification  of  the  transverse 
offset. 

The  chromatic  model  demonstrated  simulation  results  match  those  observed  by 
the  GCTEx.  The  chromatic  simulations  showed  the  reconstruction  performance  is 
degraded  over  most  of  the  spectral  range.  It  also  shifts  the  spectral  location  of  the 
best  reconstructed  image  according  to  where  the  focus  was  set.  The  chromatic  test 


case  also  showed  that  the  cone  of  missing  information  limits  reconstruction  of  high 
spatial  frequency  objects. 

The  Zernike  phase  surface  model  was  able  to  control  amounts  of  spherical,  astig¬ 
matism,  and  coma  aberrations.  The  results  showed  that  1  wave  of  aberrations  had 
little  effect  on  the  reconstruction  performance.  However,  increasing  the  amount  of 
aberrations  to  5  waves,  make  the  reconstructed  images  spatial  features  difficult  to 
distinguish.  Coma  consistently  negatively  impacted  the  reconstruction  performance 
more  than  the  other  two  aberrations  simulated.  In  designing  the  space-based  system, 
reducing  coma  should  be  held  as  priority. 

5.2  Fast  Reconstruction  Algorithm 

The  fast  reconstruction  algorithm  demonstrated  that  it  can  reconstruct  scenes  spa¬ 
tial  and  spectral  content.  The  filter  applied  by  the  reconstruction  algorithm  showed 
improved  performance  for  most  cases.  There  were  test  cases,  such  as  the  high  spatial 
frequency  object,  where  the  performance  of  the  filter  was  unacceptable.  Despite  its 
deficiencies,  the  fast  reconstruction  algorithm  is  able  to  reconstruct  strong  spatial 
and  spectral  features  but  its  performance  is  limited  due  to  the  cone  of  missing  infor¬ 
mation.  The  radiometric  performance  is  not  reliable  as  shifts  in  peaks  and  varying 
intensities  were  found.  Additionally,  the  residual  noise  artifacts  from  the  spectral 
channels  are  significant  if  left  unfiltered.  The  algorithm  would  be  well-suited  for  an 
input  to  iterative  reconstruction  algorithms. 

5.3  Transverse  Offset 

The  transverse  offset  was  determined  to  be  an  alignment  error  in  the  prism,  most 
likely  from  a  rotation  of  a  glass  about  the  normal  to  the  interface  surface  where  the 
glasses  meet.  The  total  amount  of  lateral  shift  was  measured  to  be  16.9  pixels,  cor- 


responding  to  an  alignment  error  of  20.5  arc  minutes.  Since  the  transverse  offset  is 
symmetric,  the  corresponding  correction  is  simple  to  implement  by  adding  an  ad¬ 
ditional  shift  to  the  reconstruction  algorithm  and  not  correcting  for  the  transverse 
offset  significantly  degrades  performance.  Finally,  the  adjustment  in  the  reconstruc¬ 
tion  algorithm  and  upgrade  in  electronics  marked  the  first  time  AFIT  was  successful 
in  reconstructing  an  a  measured  scene  without  manually  manipulating  the  data. 

5.4  Recommendations 

The  ground  instrument  is  primed  and  ready  for  collecting  more  laboratory  and 
held  experiments.  Simple  tests  such  as  investigating  the  effects  of  blurring  due  to  the 
prism  rotation  need  to  be  investigated.  The  effects  of  prism  rotation  blurring  will 
determine  how  fast  the  prism  can  spin  for  a  given  camera  integration  time. 

More  held  experiments  are  required  to  demonstrate  that  CTI  is  a  palatable  tech¬ 
nology.  Static  scenes  and  spectrally  interesting  scenes  are  needed  to  affirm  that  the 
instrument  is  suitable  as  a  spectrometer.  Additionally,  capturing  transient  events 
is  a  high  priority  objective  for  the  project.  The  effects  of  transient  events  on  the 
reconstruction  algorithm  need  to  be  explored. 

Prior  to  extensive  laboratory  and  held  experiments,  a  more  advanced  reconstruc¬ 
tion  algorithm  is  required  to  obtain  radiometric  precision.  However,  the  fast  recon¬ 
struction  algorithm  would  serve  as  an  excellent  “quick  look”  algorithm  that  could  be 
implemented  onto  a  digital  signal  processor  and  provide  near-real  time  reconstruction. 

Finally,  the  Zernax  and  MATLAB  combination  serves  as  an  excellent  tool  for 
future  trade  study  design  choices  for  the  spaced-based  system.  Atmospheric  effects 
or  the  slewing  of  the  space  station  can  be  incorporated  into  the  simulation  to  give  a 
representative  model  of  space-system  mission. 
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Appendix  A.  Matlab  Code 


A.l  Fast  Reconstruction  Algorithm 


%%  Min— norm  Filtered  Backpro jection  Reconstruction 
%  Chad  Su ' e 

o,  o_ 
o  o 

o, 

o 

close  all; 

%%  Load  the  cine  files  and  encoder  position 

Q. 

O 

if  ~ (exist ( 1 cameraFile ' )  &&  exist (' positionFile ') ) 

[cameraFile  cameraPath]  =  uigetf ile ( ' * . cin ' ,  'Select  the  camera  ... 
data  f ile  ' ) ; 

[positionFile  positionPath]  =  uigetfile ( ' * .mat ' ,  'Select  the  ... 
adjusted  position  file'); 

end 

if  isequal (cameraFile, 0 )  ||  isequal (positionFile, 0) 

error (' Either  camera  or  encoder  data  not  selected') 

end 

%Load  the  Phantom  Libraries  to  read  in  the  cine  files 

if  ~ libisloaded ( ' phf ile ' )  %phfile  is  the  only  library  file  really  ... 
needed  to  convert  data 

%script  was  having  troubles  locating  the  drivers  to  read  camera  data 
%so  you  will  need  to  adjust  the  directory  to  where  the  . . . 

PhantonSDK  sits 
%currDir  =  pwd; 

%cd ( ' C : \Users\Chad\Documents\MATLAB\Workspace\Phantom\PhMatlabSDK\Bin\x64\ ' ) ; 
LoadPhantomLibraries ( ) ;  %this  will  load  all  libraries 
RegisterPhantom (true ) ;  %Register  the  Phantom  dll's  ignoring  ... 

connected  cameras. 

%cd ( currDir ) ; 

end 

%%  Create  the  cine  handle  from  the  cine  file. 

%  This  script  does  not  support  a  batch  of  files 

[HRES,  cineHandle]  =  PhNewCineFromFile (fullfile (cameraPath,  cameraFile) )  ; 
if  (HRES<0) 

[message]  =  PhGetErrorMessage (  HRES  ); 

error ( [ 'Cine  handle  creation  error:  '  message] ) ; 

end 

%%  Find  the  range  and  bit  depth 

%  Get  the  first  image  and  image  range  of  the  cine  file  and  the  bit  . . . 
depth  of 

%  the  camera  data.  The  image  range  is  a  subset  of  the  total  image  count. 

%  The  image  range  is  defined  by  the  Phantom  software  to  converse  . . . 
image  file 
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%  sizes 

pFirstlm  =  libpointer ( ' int 32Ptr ’ ,  0 ) ; 
plmRange  =  libpointer (' uint32Ptr 1 , 0 )  ; 
pTotallmCount  =  libpointer ( 1 uint32Ptr ' ,  0 ) ; 

PhGetCinelnf o (cineHandle,  PhFileConst . GCI-FIRSTIMAGENO,  pFirstlm) ; 
PhGetCinelnf o (cineHandle,  PhFileConst .  GCI_IMAGECOUNT,  plmRange) ; 
PhGetCinelnf o (cineHandle,  PhFileConst . GCI_TOTALIMAGECOUNT,  . . . 

pTotallmCount) ; 
imRange  =  plmRange .Value; 
firstlm  =  pFirstlm. Value; 

lastlm  =  int 32 (double ( first Im)  +  double (imRange)  —  1); 
totallmCount  =  int32 (pTotallmCount .Value) ; 

%Get  cine  image  buffer  size 
plnfVal  =  libpointer ( 1 uint32Ptr '  , 0 ) ; 

PhGetCinelnf o (cineHandle,  PhFileConst . GCI-MAXIMGSIZE,  plnfVal) ; 
imgSizelnBytes  =  plnfVal . Value; 

%Create  the  image  range  to  be  read.  Unfortunately,  the  function  . . . 
calls  are 

%not  working  correctly  to  support  multiple  images  at  a  time 
imgRange  =  get ( libstruct ( ' taglMRANGE ' ) ) ; 
imgRange.Cnt  =  1; 

imgRange . First  =  0;  %used  to  determine  the  image  size  to  preallocate  ... 
imwrite 

[HRES,  unshiftedlm,  imgHeader]  =  PhGetCinelmage (cineHandle,  ... 
imgRange,  imgSizelnBytes) ; 

%%  Setup  the  spectral  bins,  load  the  collected  frames  then  perform  . . . 
the  2D  Fourier  Transform 

o, 

o 

N  =  50;  %number  of  spectral  bins 
%wave  =  linspace  (0 . 525,  0 . 575, N) ; 
wave  =  [0.4358335  0.546  0.576];  %Hg  Pen 

%wave  =  [0.58525  0.58819  0.59448  0.60962  0.61534  0.62173  0.62665  ... 
0.63344  ... 

%  0.64016  0.65065  0.66010  0.66780  0.67170  0.69225  0.70283  0.72452  ... 

0.74889] ;  %Ne  Pen 
%wave  =  0.546; 

%Specify  which  dataset  from  the  angular  position  make  a  full  revolution 

%imageSequence  =  int32([1163  1516  1868  2220]);%Hg  Pen 

imageSequence  =  int32 ( linspace ( 1 1 63, 255 6, N) ) ;  %Hg  Pen 

%imageSequence  =  int32 ( [234  594  950  1309]);  %Ne  Pen 

%imageSequence  =  int32 ( linspace ( 1 , 137 6, N) ) ;  %Ne  Pen 

%imageSequence  =  int32([19  109  198  287]);  %Bar  chart 

%imageSequence  =  int32 ( linspace ( 1 , 356, N) ) ;  %Bar  chart 

%Load  the  angular  position  files  and  convert  the  image  sequence  to  a 
%negative  range 

load ( fullf ile (positionPath, positionFile) ) 
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%lst  col : timestamp (MJD) ,  2nd  col: start  position,  3rd  col:  end  position 
xPixels  =  imgHeader . biWidth;  yPixels  =  imgHeader . biHeight ; 
image Sequence=image Sequence— total ImCount ; 

M  =  length ( imageSequence )  ; 

%read  the  camera  files,  transform  into  frequency  space,  and  fuse  angluar 
%data 

G  =  zeros (xPixels, yPixels,  M)  ; 
prismAngle  =  zeros (M, 1) ; 
for  m=l : M; 

imgRange . First  =  imageSequence (m) ; 

[HRES,  unshiftedlm,  imgHeader]  =  PhGetCinelmage (cineHandle,  ... 

imgRange,  imgSizelnBytes ) ; 
if  (HRES  >=  0) 

[unshiftedlm]  =  ... 

Extract ImageMatrixFromlmageBuf fer (unshiftedlm,  imgHeader) ; 
img  =  reshape (unshiftedlm ', imgHeader . biWidth,  ... 

imgHeader . biHeight ) ; 
g  =  double (img) /255; 

G ( : , : ,m)  =  fft2 (g) ; 

end 

prismAngle (m)  =  ... 

deg2rad (prismAngles (total ImCount + image Sequence (m) , 2 ) ) ;  %take  . . . 
the  first  angluar  position 

end 

%%  Setup  the  Prism  Dispersion 

o, 

o 

%load  prism  dispersion  from  Zemax  model  to  be  more  accurate 
load  ' prismDispersion . mat ' 
f 3  =  85.0;  %mm 

pixPitch  =  0.02;  %mm  per  pixel 

k  =  f3*tan (interpl (prismDispersion (:,  1)  ', prismDispersion (:, 2) ,  wave)); 

N  =  length (wave) ; 

%%  Setup  the  spatial  frequency 

o, 

o 

Fs  =  1/pixPitch;  %pixels  per  mm 
%find  the  pixel  sampling  rate 
dFx  =  Fs/xPixels; 
dFy  =  Fs/yPixels; 

%Setup  the  spatial  frequency  parameters  as  matlab  likes  it  to  properly 
%shift  the  collected/simulated  frames 
Fx  =  [ 0 : dFx : Fs/2— dFx  — Fs/2 : dFx:— dFx] ; 

Fy  =  [ 0 : dFy : Fs/2— dFy  — Fs/2 : dFy:— dFy] ' ; 

%Setup  the  Bessel  filter  for  the  reconstruction.  Use  a  meshgrid,  reshape 
%the  array,  and  for  each  spatial  frequency,  take  the  ID— FFT  wrt  lambda 
[X  Y]  =  meshgrid (Fx,  Fy)  ; 

R  =  sqrt (X. "2+Y. ~2) ; 

Mbess  =  2*pi*bessel j (0, 2*pi*R ( : ) *k*0 . 001) ; 
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Mbessel  =  reshape (Mbess, [ size (X) ,  length (wave) ]) ; 

Mfilter  =  f ft (Mbessel,  [], 3) ; 

%%  Reconstruction 

%  Ready  to  reconstruct!  For  each  spatial  frequency,  find  A  (an  MxN) 
matrix 

%  then  multiply  by  the  measured  data  cube.  Sum  all  of  the  . . . 
projections  to 

%  backpro jection .  A  1— D  filter  is  then  applied  to  acheive  the  Min— Norm 
%  solution.  After  applying  the  filter,  take  the  inverse  2D  fourier 
%  transform  to  obtain  the  hyper  cube 

R  =  1 6 . 8 66*pixPitch; 

%phase  lag  term 
%offset  =  0; 

offset  =  deg2rad (— 11 . 9242 ) ; 

%offset  =  deg2rad(45); 
whandle=waitbar (0)  ; 

Fprime  =  zeros (xPixels,  yPixels,N); 
for  n=l : N 

for  m=l:M 

phi  =  prismAngle (m) ; 

Ax  =  exp (— 2*pi*li* (k (n) *sin (phi ) +R*cos (phi+of f set ) ) *Fx)  ; 

Ay  =  exp (— 2*pi*li* (k (n) *cos (phi ) — R*sin (phi+of f set )) *Fy) ; 

%  Force  conjugate  symmetry.  Otherwise  this  frequency  ... 
component  has  no 

%  corresponding  negative  frequency  to  cancel  out  its  . . . 
imaginary  part . 

Ax (xPixels/2+1 )  =  real (Ax (xPixels/2  +  1 ) ) ; 

Ay  (yPixels/2  +  1 )  =  real (Ax (yPixels/2+1 ) ) ; 

Fprime ( : , : , n)  =  G ( : , : , m) . * (Ay*Ax) +Fprime ( : , : , n) ; 

end 

wait bar (n/N, whandle,  ' Backpro jecting .  . .  1  ) 

end 

close (whandle) ; 

%%  Apply  the  filter 

o, 

o 

mu  =  1  ; 

Fw  =  f ft (Fprime,  [  ] , 3) ; 

Fmu  =  Fw  ./ (mu  +  Mfilter) ; 

F  =  if ft (Fmu, [ ] , 3) ; 

%%  Take  the  inverse  Fourier  Transform  to  get  the  image  cube 

o_ 

o 

f  =  zeros (xPixels ,  yPixels,  N)  ; 
fprime  =  zeros (xPixels,  yPixels,  N)  ; 
for  n=l:N 

fprime ( : , : , n)  =  if ft2 (Fprime ( : , : , n) ) /M; 
f (  :  ,  :  ,  n) =real (if ft2  (F  (  :  ,  :  ,  n) ) ) /M; 

end 
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A. 2  Mercury  Pen  Lamp  Simulation 


%%  Mercury  (Hg)  Pen  Lamp  Simulation  of  GCTEx 
%  Capt  Chad  Su'e 

o,  o. 
o  o 

o, 

o 

clear  all 
close  all 

%%  Load  and  define  simulation  parameters 

%  Load  common  simulation  parameters  from  simParams  class  and  define 
%  specific  parameters  for  this  simulation 

%Load  default  simulation  parameters  from  class 
hgSim  =  simParams (); 

%%  Load  &  Initialize  Zemax 

o. 

o 

zDDEInitO;  %need  an  error  check  to  determine  if  Zemax  is  open 
zLoadFile (hgSim . modelPath)  ; 
zPushLens ( 1 ) ; 

%%  Load  Image  Simulation  Parameters 

o, 

o 

%  TODO:  implement  a  shape  vs  time  profile 
imgSim  =  imgSimParams ( ) ; 

%imgSim. ISM.INPUTFILE  =  ' ALPHABET . IMA '  ; 

%imgSim.  ISM_FIELDHEIGHT  =  0.4; 
imgSim.  ISM.INPUTFILE  =  ' CIRCLE . IMA ' ; 
imgSim.  ISMJFIELDHEIGHT  =  0.02; 
imgSim . updateSettings ()  ; 

%%  Set  up  prism  angles 

o. 

o 

%For  reconstruction,  were  are  going  to  use  25  wavelength  bins,  so  ... 
break  the 

%prism  rotations  into  25  equal  rotations 

%TODO:  integrate  the  discrete  nature  of  the  encoder 
%bins  =  100; 
bins  =  36; 

prismSteps  =  360/bins; 

prismAngles  =  linspace (0, 360— prismSteps, bins) ; 
hgSim . numAngles  =  length (prismAngles ) ; 


%%  Configure  the  wavelength  profile 

o, 

o 
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wave  =  [0.4046565  0.4358335  0.5460750  0.5769610  0.5790670];  %um 
wavelntensity  =  [4400  10000  10000  1100  1200]; 

wavelntensity  =  wavelntensity/max (wavelntensity ) ; 

hgSim . numWave  =  length (wave ) ;  %keep  track  in  the  simulation  class 
%%  Start  the  simulation 

o, 

o 

%The  text  file  is  a  large  file  size,  so  rather  parsing  a  text  file,  . . . 

Zemax 

%allows  for  the  simulation  results  to  be  output  the  sim  as  a  rgb  bitmap. 

%The  sim  settings  will  always  output  a  each  wavelength  iteration  . . . 
onto  the 

%blue  channel  and  the  text  file  is  completely  ignored.  For  each  ... 
wavelength 

%iteration,  the  sim  output  is  normalized  and  then  scaled.  For  each  angle 
%iteration,  the  wavelengths  results  are  then  summed  and  saved. 

%  simResults  holds  the  sim  outputs  for  each  wavelength  for  every  angle 
g  =  zeros (hgSim . xPixels ,  hgSim. yPixels,  hgSim . numAngles ) ; 
simResults  =  zeros (hgSim . xPixels ,  hgSim . yPixels ) ; 
tempResult  =  zeros (hgSim . xPixels ,  hgSim . yPixels ) ; 

whandle=waitbar (0)  ; 
for  j  =  1 : hgSim . numAngles 
%set  the  prism  angle 
zSetMulticon (1,1, prismAngles (  j  )  )  ; 
for  k  =  1 : hgSim. numWave 

waveData  =  zSetWave ( 1 , wave (k) , 1 )  ; 

status  =  zGetTextFile ( [hgSim. simTempResultPath  ... 

'\tempSimResult.txt'],  'Sim',  hgSim. settingsFile,  ... 
hgSim . flag) ; 

tempResult=imread ( [hgSim. simTempResultPath  . . . 

' \tempSimResult . bmp ' ] ) ; 

%get  sim  results  from  blue  channel,  normalize  and  then  scale. 
simResults=simResults+double (tempResult ( : , : , 3) ) /255*wavelntensity (k) ; 

end 

g(:,:,  j)  =  simResults; 

simResults  =  zeros (hgSim . xPixels ,  hgSim. yPixels) ; 
waitbar ( j /hgSim. numAngles , whandle, ' Simulating ...'); 

end 

close (whandle) ; 

%%  Save  simulation  output  and  close  DDE  server 

o, 

o 

prismAngles  =  deg2rad (prismAngles) ; 

save ( [hgSim. simResultsPath  ' \simout .mat ' ] , ' g ' , ' prismAngles ' , '—mat ' ) ; 
zDDEClose ()  ; 
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A. 3  Simulation  Parameters 


classdef  simParams 

%  %  define  the  default  properties  of  the  simulation 
properties (GetAccess  =  'public',  SetAccess  =  'public') 

%%Directory  &  Filename  Path  Parameters 
modelPath  =  ... 

' C : \Users\Chad\Documents\MATLAB\Workspace\CTEx\Simulation\Models\GCTEx_Bosti 
simResultsPath  =  ... 

' C : \Users\Chad\Documents\MATLAB\Workspace\CTEx\Simulation\hgSimResults ' ; 
simTempResultPath  =  'C:\Users\Chad\Documents\ZEMAX\lMAFiles'; 
settingsFile  =  ... 

' C:\Users\Chad\Documents\ZEMAX\Configs\MySIM.CFG'  ; 

%%Prism  Parameters 

prismSpeed  =  15;  %revolutions/second 
numAngles ; 

%%Wavelength  Parameters 
numWave ; 

%%Image  Simulation  &  Camera  Parameters 

INPUTFILE  =  ' CIRCLE . IMA ' ;  %  The  input  file  in  which  to  simulate 

xPixels  =  512; 

yPixels  =  512; 

pixPitch  =  0.02;  %units:  mm 

flag  =  1;  %img  sim  settings:  0— default,  1— settings  file,  ... 

2— prompt 

%%Encoder  Parameters  %maybe  add  prism  angles 
encoderBitResolution  =  12; 
encode rRe solution; 


end 

properties (Constant  =  true) 
%add  any  constants  here 


end 

methods 

%  methods,  including  the  constructor  are  defined  in  this  block 
function  obj  =  simParams () 

%  class  constructor 

ob j . encoderResolution  =  2 "obj . encoderBitResolution; 

end 

%Define  what  Zemax  model  to  load 
function  obj  =  setZemaxModel (ob j , f ilePath) 
if (nargin  >  0 ) 
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ob j . modelPath  =  filePath; 

else 

%add  a  filepath  chooser 

end 

end 

%Define  the  folder  directory  to  store  simulation  results 
function  obj  =  setResult sDir (ob j , dir ) 
if (nargin  >  0 ) 

ob j . simResultsPath  =  dir; 

else 

end 

end 

end 

end 


98 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


A. 4  Image  Simulation  Parameters 


classdef  imgSimParams 

%%Define  the  default  properties  of  the  simulation 
properties (GetAccess  =  'public',  SetAccess  =  'public') 

%  ISM_INPUTFILE :  The  input  file  name. 

%  ISM_FIELDHEIGHT :  The  Y  field  height. 

%  ISM_OVERSAMPLING :  Oversample  value.  Use  0  for  none,  1  for  2X,  2  . . . 
for  4x, etc . 

%  ISM_GUARDBAND :  Guard  band  value.  Use  0  for  none,  1  for  2X,  2  for  . . 
4x,  etc. 

%  ISM_FLIP:  Flip  Source.  Use  0  for  none,  1  for  TB,  2  for  LR,  3  for  . . 
TB&LR . 

%  ISM_ROTATE:  Rotate  Source:  Use  0  for  none,  1  for  90,  2  for  180,  3  . 
for  270. 

%  ISM_WAVE:  Wavelength.  Use  0  for  RGB,  1  for  1+2+3,  2  for  wave  #1,  3 
for  wave#2,  etc. 

%  ISM.FIELD:  Field  number. 

%  ISM_PSAMP:  Pupil  Sampling.  Use  1  for  32x32,  2  for  64x64,  etc. 

%  ISM_ISAMP:  Image  Sampling.  Use  1  for  32x32,  2  for  64x64,  etc. 

%  ISRLPSFX,  ISFLPSFY:  The  number  of  PSF  grid  points. 

%  ISM_ABERRATIONS :  Use  0  for  none,  1  for  geometric,  2  for  diffraction 
%  ISM_POLARIZATION :  Use  0  for  no,  1  for  yes. 

%  ISM_SHOWAS:  Use  0  for  Simulated  Image,  1  for  Source  Bitmap,  and  2  . 
for  PSF  Grid. 

%  ISM_REFERENCE :  Use  0  for  chief  ray,  1  for  vertex,  2  for  primary  . . . 
chief  ray. 

%  ISM.SUPPRESS :  Use  0  for  no,  1  for  yes. 

%  ISM_PIXELSIZE :  Use  0  for  default  or  the  size  in  lens  units. 

%  ISM_XSIZE,  ISM_YSIZE:  Use  0  for  default  or  the  number  of  pixels. 

%  ISM_OUTPUTFILE :  The  output  file  name  or  empty  string  for  no  output 
f  ile  . 

settingsFileName  =  ... 

' C:\Users\Chad\Documents\ZEMAX\Configs\MySIM.CFG'  ; 
ISM.INPUTFILE  =  ' CIRCLE . IMA '  ; 

ISM.FIELDHEIGHT  =  0.02; 

ISM_OVERSAMPLING  =  uintl6(0); 

ISMJGAURDBAND  =  uintl6(0); 

ISM.FLIP  =  uintl 6  ( 0 ) ; 

ISM_ROTATE  =  uintl6(0); 

ISM.WAVE  =  uintl6 (2) 

ISMLFIELD  =  uintl6 (1) 

ISMLPSAMP  =  uintl6 (1) 

I SM_I SAMP  =  uintl6 (1) 

I SM_P  SFX  =  uintl 6 (3) 

I SM_P  SFY  =  uintl 6 (3) 

ISM_ABERRATIONS  =  uintl6(2) 

ISMLPOLARIZATION  =  uintl6(0) 

ISM.SHOWAS  =  uintl 6 ( 0 ) 

I  SM_REFERENCE  =  uintl6(l) 
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ISM.SUPPRESS  =  uintl6 (0) 

I  SM_P IXELS I ZE  =  0.02 
I SM_XS  I  ZE  =  uintl6 (512) 

I SM_YS I ZE  =  uintl6 (512) 

ISM.OUTPUTFILE  =  'tempSimResult.bmp' 

end 

methods 

%  methods,  including  the  constructor  are  defined  in  this  block 
function  obj  =  imgSimParams ( ) 

%  class  constructor 


end 

function  obj  =  updateSettings (ob j ) 

zSetSettingsFile (obj . settingsFileName, 
obj . I SM_INPUTFILE ) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISMJFIELDHEIGHT) ; 
zSetSettingsFile (obj . settingsFileName, 
obj . I SM_OVERS AMP LING) ; 
zSetSettingsFile (obj . settingsFileName, 
obj  .  ISMJ3AURDBAND)  ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISM.FLIP) ; 

zSetSettingsFile (obj . settingsFileName, 
obj  .  ISM_ROTATE) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . I SM_WAVE ) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISM.FIELD) ; 

zSetSettingsFile (obj . settingsFileName, 
obj  .  I  SM_P  SAMP  )  ; 

zSetSettingsFile (obj . settingsFileName, 
zSetSettingsFile (obj . settingsFileName, 
zSetSettingsFile (obj . settingsFileName, 
zSetSettingsFile (obj . settingsFileName, 
zSetSettingsFile (obj . settingsFileName, 
zSetSettingsFile (obj . settingsFileName, 
obj  .  I SM_SH0WAS ) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISM.REFERENCE) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISMLSUPPRESS) ; 

zSetSettingsFile (obj . settingsFileName, 
obj  .  I  SM_P  IXELS  I  ZE )  ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISM.XSIZE) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . I SM_Y  SIZE) ; 

zSetSettingsFile (obj . settingsFileName, 
obj . ISM.OUTPUTFILE) ; 


' ISM.INPUTFILE  '  ,  ... 

' ISM_FIELDHEIGHT  '  ,  ... 

' ISM.OVERSAMPLING ' ,  ... 

'  ISM_GAURDBAND ' ,  ... 

'  I SM_FLIP  '  ,  ... 

' I SM_R0TATE ' ,  ... 

'  ISMJWAVE  '  ,  ... 

'  ISM.FIELD  '  ,  ... 

'  I  SM_PSAMP  '  ,  ... 

'  I  SM_I  SAMP  '  ,  obj  .  I  SM_I  SAMP  )  ; 

'  I  SM_P  SFX  '  ,  ob  j  .  ISM.PSFX)  ; 

'  I  SM_P  SFY  '  ,  ob  j  .  I  SM_P  SFY )  ; 

' ISM.ABERRATIONS ' , obj  .  ISM_ABERRATIONS ' 
’ I S NLP 0 LARI ZATION ’ , obj  .  ISM_POLARIZATI( 
' I SM_SH0WAS ' ,  ... 

'  ISM_REFERENCE ' ,  ... 

' ISMLSUPPRESS ' ,  ... 

'  I  SM_P  IXELS  I  ZE  '  ,  ... 

' I SM_XS I ZE ' ,  ... 

' I SM_YS I ZE ' ,  ... 

'  ISM.OUTPUTFILE ' ,  ... 


end 


end 


end 
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Appendix  B.  GCTEx  Electrical  System  Block  Diagram 
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